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I. INTRODUCTION 



The understanding of hadronic systems is at the heart of numerous theoretical studies these last 20 years. 

While the description of experimental processes at very high energy can rely on perturbation theory, the study 
of bound states necessitates the use of non-pcrturbativc methods. Concerning the latters, the nonperturbative 
renormalization procedure and the invariance of the system under gauge transformations are the properties 
of particular interest. Among the various frameworks used over the past to tackle these problems, we shall 
concentrate on Light Front Dynamics (LFD) [1] which exists in two field-theoretical forms: standard LFD [2] 
and explicitly Covariant Light Front Dynamics (CLFD) [3]. Standard LFD deals with the state vector defined 
on the plane t + z = 0, while in CLFD this plane is given by the invariant equation iv-x = 0, where u) is an 
arbitrary four- vector with = 0. The particular choice ui = (1, 0, 0, —1) turns CLFD into standard LFD. We 
shall sec below that the explicit covariance of CLFD is a very powerfuU tool to address the above questions in 
a transparent way. 

Following the pioneering work of Wilson [4], the perturbative renormalization in QED has been studied 

both in LFD [5-7] and CLFD [8]. In these papers the fcrmion self-energy is investigated, and the necessary 
counterterms needed to renormalize it are found. To generalize the renormalization technique for nonpertur- 
bative processes, it is necessary to consider the same questions in terms of the Fock decomposition of the state 
vector. This has been already done for scalar models in LFD [9-11] and CLFD [12]. Such toy models are 
important to settle general equations and test numerical procedures. At the same time they are too simple 
because of, firstly, their super-renormalizability, and, secondly, the absense of contact interactions appearing for 
non-zero spin particles in LFD and CLFD. The next step therefore is an investigation of fermionic composite 
states [9,13-16]. Two important properties have emerged from these calculations. 

i) The possible Fock sector dependence of the (mass) countcrtcrm, which has been discussed in Ref [9] 
for LFD. The necessity to remove, by hands, the counterterm in the equation which defines the last Fock state 
component considered in the calculation is now well understood [11,12]. This (infinite) counterterm should be 
balanced by a (also infinite) contribution arising from higher Fock state components not taken into account due 
to the truncation of the Fock state expansion. 

ii) The appearance of non-local counterterms needed to recover known symmetries and to renormalize 
the theory [6]. As has already been shown in CLFD [8], the knowledge of the structure of the counterterms as 
functions of the light front plane orientation is necessary to understand how non-local counterterms appear in 
LFD. We shall come back extensively to this question in the following. 

We consider in the present study a fermion F coupled to bosons b with spin either zero or one. We 
restrict ourselves to the approximation where only the first two Fock sectors (j-F) + \Fh)) in the state vector are 
retained. Such a model is instructive in many aspects: [i) since we arc dealing with fermions, we have to take 
into account contact interactions. This, also, implies non-trivial contributions arising from the counterterms 
needed to renormalize the theory; [ii) we can investigate the case of gauge theories; {Hi) the procedure can 
be checked by comparison with results of previous calculations based on perturbative methods, as mentioned 
above. It can thus be extended with some confidence to calculations involving higher Fock state components. 

In a first step, we suppose the boson h to be spinless, that just corresponds to the well-known Yukawa 
model. This gives us an opportunity to develop our method and show the peculiarities of CLFD. We reproduce 
the results obtained in this model in Ref. [14] in standard LFD. 

We then study the system where the constituent boson is a gauge particle. In order to investigate the 
infiuence of the gauge on the state vector, we consider two difi'erent gauges: the Feynman and the Light-Cone 
(LC) ones. We will see that the number of spin components of the two-body wave function depends on the 
gauge. In the Feynman gauge the wave function is determined, in general, by eight components, whereas in 
the LC gauge it contains four components. However, in the simple model developed in the present paper 
where we restrict ourselves to the \F) + \Fb) Tamm-Dancoff truncation, three non-zero components remain 
for the Feynman gauge and two ones for the LC gauge. Moreover, after imposing a certain condition on the 
counterterms, only one non-zero component of the two-body wave function survives for each gauge. 

The plan of the paper is as follows. In Sec. II we present our general framework. In Sec;. Ill the need for 
a specific extra counterterm in LFD is justified. In Sec. IV we calculate the state vector and the counterterms 
for the scalar boson case. In Sec. V we solve the problem for the gauge boson in the Feynman gauge. The same 
problem, but for the LC gauge, is solved in Sec. VI. The results arc summarized in Sec. VII, where we discuss 
also the questions of gauge invariance and comparison with perturbative calculations. Sec. VIII contains our 
concluding remarks. Technical derivations are rounded up in Appendices. 
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II. GENERAL CONSIDERATIONS 



A. Exact Lagrangian for scalar bosons 

In order to clearly settle the starting point of our study, we shall recall in this section some known results. 
For simplicity, we detail only the scalar boson case. The case of vector bosons can be deduced very easily. 
We shall start with the bare Lagrangian written as 

£ = *o [i^ - mo] *o + ^ [d^^od^^o - Mo^o] + 5o*o*o$o, (1) 

where ^ = j-d = ^"di,, (^o) is the bare fermion (boson) field, nio (/io) is the bare fermion (boson) mass, 
and go is the bare coupling constant. Going over from bare to dressed quantities by means of the standard 
substitutions 

*0 = ^2^'*, $0 = ^3^'$' 5^1=50^2^3^', (2) 

where g is the physically observed value of the coupling constant, and introducing also the physical masses m 
and fj,, we get 

£ = §[i,0t-TO]* + i [d^<^d''^ - i^^<P^] + 5*** 

+ {Z2 - 1)* [i^ - m] * + i^^-lll [d^^d''^ - ^^$2] 

+ (5m** + + g{Zi - 1)**$. (3) 

From Eqs. (l)-(3) it is easy to find that 

5m = Z2{m-mo), J/i^ = .^3(a*^ - Mo)- (4) 

The countcrtcrms Zis, Sm, and (J/i^ can be foimd from the definitions of the observed quantities g, 
m, and /i^. Namely, dm and Sfi^ are defined by the positions of the one- fermion and one-boson propagator 
poles. The additional requirement that the residues of these propagators at their poles equal i allows to fix the 
constants Z2 and Z3. Finally, the constant Zi is determined from the condition of the coincidence of the on-shell 
dressed fermion- fermion-boson vertex at zero boson four-momentum with the physical coupling constant. 

B. Perturbative vs Fock state decomposition 

As we already mentioned in the Introduction, a natural framework to study composite systems in LFD 
is the Fock state decomposition. Renormalization procedures, on the other hand, are well under control in 
perturbative Quantum Field Theory using the 4-Dimcnsional Fcynman (4DF) diagrammatic expansion. In the 
latter approach, the composite fermion mass is defined as the pole of the fermion 2-Point Green's Function 
(2PGF) which is obtained by summing up the chain series of the iterations of the fermion self-energy (a sum of 
all irreducible Feynman diagrams calculated up to a given order of coupling constant), S4£)ir(^, mo), according 
to 



it % % 

12 + [-i'^iOFi'^, mo)] + . . . = = -rz T 

p — mo p — mo p — rriQ p — mo — i-,4DF{p,'mo) 



(5) 



In LFD state vectors are defined as cigcnstates of the corresponding light-front Hamiltonians. By using 
the Fock decomposition of the state vector, the eigenstate problem can be easily reduced to that of solving a 
system of linear homogeneous equations which involve various Fock components. Note that the Fock decom- 
position is made in terms of the number of Fock components rather than the coupling constant powers. Each 

2 

Fock component effectively brings irreducible contributions of order a or higher, with a = f;;, to the fermion 
propagator. The summation of these contributions to all orders in a is done implicitly when solving the system 
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of linear equations mentioned above. The mass of the composite particle is found from the condition that this 
system has a nontrivial solution. Such a procedure is an analogue, in LFD, of the summation of the chain 
series for the 2PGF in 4DF approach. Below wc shall derive and solve this system of equations, for the Fock 
decomposition restricted to the components \F) + \Fb), for scalar and gauge bosons in both the Feynman and 
LC gauges. 

C. Renormalization scheme 

Before going into more detailed and technical considerations, it is necessary to clarify in which scheme 
we shall work to renormalize the system under study. A nice exposition of the standard derivation can be found 
in [17,18]. The following two schemes are usually considered in perturbative calculations: 

• Bare renormalization scheme. The starting point is the bare Lagrangian with bare masses and coupling 
constants, as given by Eq. (1). These masses and coupling constants are determined from physical condi- 
tions. For instance, the bare fermion mass is chosen so to get a pole of the 2PGF at the physical fermion 
mass. The bare mass expressed in terms of the physical one is in general infinite and must be regularized. 

• Dressed renormMization schem,e. The starting point here is a rcnormalizcd Lagrangian (3) defined in 
terms of (finite) physical masses and coupling constants, and including the counterterms which are thus 
infinite and must be regularized. 

The two regularization schemes should of course give the same physical observables provided no approximations 
have been made. The latter never happens in practice: one cither makes a perturbative expansion, or truncate 
Fock space in LFD, or both. The equivalence between the two schemes can be directly checked in the lowest 
order of perturbation theory, as explained in any text book on Quantum Field Theory. In LFD, this question 
was partially adressed in Refs. [4,14]. We have also tackled it within CLFD using perturbative expansions [8]. 

Nonperturbative calculations however hardly admit such a check. In the spirit of the Fock state decom- 
position, one should make sure that the one-particle states arc as close as possible to the physical ones. In 
our case this implies that the basic one- fermion state \F) and the full state \F) + \Fb) correspond to the same 
physical mass m, and we thus should consider the mass counterterm 5m. 

The counterterm proportional to Z2 — 1 in Erj. (3) renormalizes the fermion field. However, since we arc 
working here with the state vector rather than with particle fields, it is more appropriate to deal directly with 
the state vector normalization factor A'' instead of Z2. For this reason, we will exclude the counterterm with Z2 
from our consideration. The normalization factor N does not enter the eigenstate equation for the state vector 
and can be calculated after the solution is found. 

Our task is simplified for the Fock decomposition truncated to the components |F) -|- \ Fb). In this case, 
the boson state vector is still free, therefore, Z3 = 1 and J/i^ = 0. 

Finally, the counterterm with Zi is related to the coupling constant renormalization. Since in the present 
paper we do not calculate physical observables like electromagnetic form factors or scattering amplitudes, we 
may deal with the unrenormalized coupling constant and put Z\ = \. For convenience, we shall call this coupling 
constant g. 

So, as far as the counterterms in the Lagrangian are concerned [the second and the third lines in Eq. (3)], 
it is enough for solving the nonperturbative problem in the two-body approximation, to retain explicitly the 
term with 5m only, while those containing Z\, Z2, Zs, and can be dropped out. The quantity dm defined 
now as 

Sm = m — mo, (6) 
that is a shift of mass caused by the interaction. 

III. THE NEED FOR A NEW COUNTERTERM IN LFD 

As we shall see explicitly in the following applications, it turns out that one has to introduce in LFD 
another counterterm which is not reduced to any of those considered above. 
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To explain its role, let us calculate the light-front 2PGF in the dressed renormalization scheme, using 
the two-body approximation. Below in this section we obtain a formal expression for the 2PGF, then find the 
pcrturbativc (up to terms of order g^) fcrmion-boson on-shcll scattering amplitude M^^^ and demonstrate that 
without an additional counterterm such an amplitude would be dependent on the light front orientation (i. e. 
on the four- vector w). The dependence of approximate physical CLFD amplitudes on ui is equivalent to the 
rotational symmetry violation in standard LFD. 

Since we are not dealing here with an eigenstate problem, we should retain the constant Z2 in the 
Lagrangian (3), while putting Zi = Zs = 1, Sfi"^ = 0, as explained above. So, the Lagrangian we consider in 
this section is 

£ = ^[z^-m]*-|-^ [d^^d"^ - -I- 5**$ + {Z2 - 1)§ [i ^ - m] * ^m**. (7) 

In the following we will use the fact that perturbative expansions of the quantities Z2 — 1 and 6m start with 
terms of order 9^. 

Denote by pi (p'l) the initial (final) fcrmion four- momentum; the corresponding boson four-momenta 
are ki {k{). In the lowest (second) order in g the on-energy-shell scattering amplitude M^^) is defined by the 
sum of the two CLFD diagrams shown in Fig. 1. The rules of the CLFD diagram technique are exposed in 
detail in Ref. [3]. As usual, straight solid, wavy, and dashed lines correspond to fermions, bosons, and spurious, 
respectively. The internal fermion line with a large full blob in the second diagram is the so-called contact 
term^. The amplitude M^^^ was calculated in Ref. [3]: 



M(2) = g^u{p[) 



p2 2{lu-p) 



(8) 



where p = pi + ki = p'l + k'l, r = {m^ — p'^)/2{uj-p). The two addcndums in the square brackets on the r. h. s. 
of Eq. (8) just correspond to the contributions from the two diagrams shown in Fig. 1. It is easy to see that 
the w-dependent terms in the sum cancel each other, the whole amplitude M^^^ being w-independent: 

M(2) = ig^u{p[)D{p)u{pi), (9) 

where 

Dip)=i4^ (10) 
p'' — 

is the 2PGF in the lowest (zeroth) order of perturbation theory. Note that it has the same form as in the 4DF 
approach, as it should. 

The exact light-front 2PGF can be obtained by formal summing up the perturbation series. Such a 
procedure is quite similar to that used in standard QED and leading to the Dyson equation. The difference is 
that LFD calculations involve contact terms mentioned above. 

Let us define the mass operator ~M(j}) (the sign minus is introduced for convenience) which is a sum of 
all irreducible light-front diagrams containing no one-fermion intermediate states or, in other words, including 
only intermediate states with bosons. The equation for the exact 2PGF I?(p) is obtained by infinite iterations 
of the kernel being a sum of the three terms: the vertices {Z2 — 1)(^ — m), dm, and the mass operator —A4{p). 
This equation reads 



V{p) = D{p) + V{p)[i{Z2 - l){l^-m) + iSm - iM{p)]D{p). (11) 

il solutioi 

V{p)=D{p)-_ 



It has the formal solution 

1 



1 - i[{Z2 - l){fi-m) + 6m- M{p)]D{p) 



= i (121 

p'^-m'^ + {Z2-l){p^-m'^) + [6m-M{p)]{fl + m)' ^ ' 



^In Ref. [3] the eoritact terms were represented by lines with crosses. Here we changed this notation, reserving the cross 
symbol for counterterms (see below). 
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In the g'^ (or two-body) approximation the mass operator is given by the sum of the two diagrams 
shown in Fig. 2. The diagrams (a) and (b) taken with the minus sign each are the light-front fermion self- 
energy, — S(p), and a contact term on the fermion Hnc, — E^c(p), respectively. Note that the amplitudes of 
these diagrams include integrals over bosonic momentum. Because of different phase space volumes, it is not 
judicious to combine the contact term with the usual light-front fermion propagator in the self-energy diagram 
into the unique term D{p), like it has been done in Eqs. (8) (10). So, the diagrams (a) and (b) in Fig. 2 
must be treated separately. The mass operator M {p) is a light-front analogue of the 4DF fermion self-energy 
'^ADwi'^ifn) calculated in the dressed renormalization scheme. 

Now let us expand Eq. (12) in powers of g up to terms of order g"^. We have 



V{p) « D(^)(p) . z4±^ - ^4±^^ f(Z, - m -m)+Sm- M^^\p)] 



(13) 



with 



Applying the CLFD graph technique rules [3] to the diagrams shown in Fig. 2, we get 

S(p) = --f-^ / {fi + m)9{uj-q)6{q'^ - m^)d'^q e{uj-k)5{k'^ - ii^)d'^k6^'^\p + lot - k - q)- 

J '. 



dr 



(14) 



iO 



{2nf 



m) 9[Loip - k)] d^k _ 



2u-{p — k)T 



^(p2)+B(p2)i+C(p2) 
m 



m i/b 



uj-p 



and 



where 



S/c(p) 



1 



(27r)3 j 2w(p - k) 2ek 



d^k 0^ m il) 
— = g^Cf, 



u)-p 



T 



2ujip-k) 



(15) 



(16) 



(17) 



and Ek = ^/k^ + Ji^. The functions A{p'^), B{p'^), C(p^), and the constant C/c are given in Appendix Bl. 
Note that the self-energy T,{p) depends on the position of the light front [8], since it contains the contribution 
proportional to j^. The same contribution was found earlier in Ref. [6]. In Eq. (B9) from Ref. [6] it looks as 

that is just the value of ^ for u> = (1, 0, 0,-1). Substituting the decompositions of E(p) and S/c(p) in 
terms of A{p^), B{p^), C{p'^), and Cfc, as defined by the r. h. s. of Eqs. (15) and (16), into Eq. (14) we can 
write 



M^'\p) = g' {a{p') + B{p')l + [C(/) + Cf^] ^1 . 

1^ m ui-p ) 



(18) 



The fermion-boson scattering amplitude calculated in the two-body approximation discussed here can 
be obtained by sandwiching the 2PGF (12) between the bispinors u{Pi) and u{pi). Up to terms of order 5^ 
inclusive it is 



M(4) =ig^u{p[)V^^\p)u{pi). 
Substituting Eq. (18) into Eq. (13) and then into Eq. (19), we find after some algebra: 



(19) 



(20) 



where 
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^ [Z,-2- g^g(p^)/m](^ + m) + - g^[.4b^) + Bjp')] \ 

— m? / ' 

MW=,^^^^1±%.(K)^«(pO- (22) 

The constants Z2 and are found in a standard way (we will specify the corresponding procedure below) and 
do not include any dependence on the orientation of the light front uj. As a result, the quantity is also 

a;-indcpcndcnt. On the contrary, explicitly depends on u) through the term u{pi) i/ju{pi). So, the whole 

amplitude M^^^ depends on w as well, that evidently calls for a certain additional counterterm (except for the 
previously introduced Z2 and 5m) to eliminate this nonphysical w-dependence. 

On can also come to the same conclusion from the requirement that the 2PGF (13) has a simple pole at 
^ = m with the residue equal to i. In other words, there must be 

^ W , =-, + --- = i^ 2+---' 23 

p^rn p — m p'^ — m^ 

where the dots designate finite terms. Comparing Eqs. (23) and (13), we have to demand that 

li^ \ {l^+m)[{Z2-m-m) + 5m-M^^\p)]{^ + m) \_^_ 
p2^m^ \ p^ -m^ J ' 

The expression in the curled brackets in Eq. (24) is a matrix. The condition (24) is imposed on any of its matrix 
element. Substituting Eq. (18) into Eq. (24) and using the Taylor expansions A{p'^) = A + A' ■ [p^ — m^) + . . ., 
Bip^) = B + B' ■ {jp^ - m^) + C{p'^) = C -\- C ■ {p^ - m^) + ... near the point p^ = m^, we can rewrite 
Eq. (24) in the form 

o (j.^ ^ 5m-g\A + B + C + Cfc) 



P'^ — 

+ [m{Z2 - 1) - g'^B - 2g^m^{A' + B' + C) + 6m-g'^{A + B) 



+ [m{Z2 - 1) - g'^B - 2/m2(A' + B' + C")]- + g\C + QJ— = 0. (25) 

m uj-p 



The absense of the pole term in Eq. (25) leads to 



.2/ . , D , , ^ ^ u{p)M'^'^>{p)u{p) 



5m = g\A + B + C + C}c) 



2m 



(26) 



that coincides with the standard perturbative expression. The last equality immediately follows from Eq. (18) 
and the identities u{p)u{p) = 2m and u{p)j^u{p) = p'^u{p)u{p)/m valid at p'^ = m^ . However, it is impossible 
to choose Z2 in Eq. (25) so that all the other coefficients at the independent matrices 1, ^, and ^ would be 
zero, unless C + C/c = 0. As is seen from Eqs. (B4c) and (B4d) in Appendix B 1, the sum C + C/c, generally 
speaking, differs from zero. Hence, Eq. (25) can not be satisfied. This again indicates the need for another 
counterterm in the Lagrangian (7). 

It is easy to see that the counterterm required must have the following structure (in momentum repre- 
sentation): 

^ - 1) , (27) 

where Z^ is a constant to be determined. We have supplied this constant with the subscript "w" in order to 

point out its origin, but it does not in fact depend on lo. The precise form of Eq. (27) is chosen so that this 
counterterm does not contribute when sandwiched between two bispinors of equal momenta. 
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After introducing the new counterterm, the quantity 6m — M.^^\p) in Eqs. (13) and (24) must be replaced 

by 



5m-A^(2)(p) + Z^(^-l 



that is equivalent to the following two substitutions: 

5m^Sm-Z^, g^[C{p^) + Cjc] ^ g^[C{p^) + Cfc] ^ Z^, (28) 
made in all equations obtained above in this section. Now, instead of Eq. (25) we get 

Sm-g^{A + B + C + Cfc) 



2m(^ + m) 2 

P ■■- — 

+ [m{Z2 - 1) - g^B - 1g^m^{A! + B' + C) + 5m - Z^ - g'^{A + B)] 

+ [m{Z2 - 1) - g^B - 2g^m\A' + B' + C')]- + [g\C + Cf^) - Z^]^ = 0. (29) 

m u)-p 



Prom here we reproduce for Sm the same formula (26), while for Z2 and Z^ we find 

Z2 = 1 + ^B + 2mg\A' + B' + C')=l + ^{ u{p)^^-^u{p) 
m 2m I op 



(30) 



Z^=g^{C + Cf,). (31) 

In differentiating, A4^'^''{p) is implied to depend on both explicitly and through p"^ = /p- /p- Note also 
that Z2 has the same form as in the standard renormalization technique. So, after introducing the additional 
counterterm (27) the behaviour of the fermion propagator near its pole is indeed governed by Eq. (23). 

One can easily check that the counterterm (27) eliminates the w-dependence of the scattering ampli- 
tude (20)-(22). This nonphysical w-dependence is concentrated in the term M.'^f^ given by Eq. (22). The 
substitutions (28) with Z^ taken from Eq. (31) transform Eq. (22) into 

p^ —m^ uj-p p^ — m^ (jo-p 

since, by definition, C = C{m^). As proved in Appendix Bl, C{p^) does not depend on p^ [in contrast to 
the functions A{p'^) and B{p'^)]. Therefore M^^ = and the whole scattering amplitude (20) turns out to be 
w-independent, as it should. 

We have constructed the additional counterterm (27) from perturbative considerations. However, as will 
be shown below, it retains its form for the nonpcrturbativc problem discussed in the present paper. Moreover, 
its structure is the same for both the scalar and gauge boson cases. In the next section and in Appendix A we 
will explain how to introduce it in the light-front Hamiltonian. 

IV. SCALAR BOSON 

In this section wc give a solution for the state vector of a system composed of a spin 1/2 fermion coupled 
to a scalar boson, in the Fock decomposition restricted to the \F) + \Fb} approximation. A general derivation 
of the CLFD fermion-boson Hamiltonian for the scalar boson case is presented in Appendix A 1. 

A. Contact terms 

As already mentioned in Sec. II, a peculiarity of LFD is the appearcncc, for particles with spin, of the 
so-called contact interactions [2,3]. These interactions arise from the elimination of two non-dynamical degrees 
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of freedom, for the fcrniion field, as recalled in Appendix A. In the case of spinor-scalar interactions, with the 
initial interaction Lagrangian of the form Z/'"*(a;) = gxptpip, where xp and ip arc, respectively, the free fermion 
and boson fields, this contact interaction generates a FFbb cx vertex in addition to the usual FFb vertex. In 
the presence of counterterms needed to renormalize the theory, like the standard mass counterterm oc Sm'tptp, 
the elimination of the non-dynamical degrees of freedom generates in addition two new contributions to the 
light-front Haniiltonian. The first one proportional to Sm/^ originates from two self-interaction vertices ^, the 
second one proportional to g6m appears from the product of the self-interaction times the ordinary FFb vertex. 
Taking into account the specific counterterm discussed in Sec. Ill and introducing, for convenience, a shifted 
mass counterterm Sm' defined by 

5m' = 6m-Z^, (33) 

we obtain the final form of the light-front interaction Hamiltonian. The details of the derivation are given in 
Appendix A 1. We thus get 

fl''"*(a;) = Hi{x) + Hu{x) + H^ix) + H2c{x) + H^^{x) + Hsix) (34) 



with 



Hi{x) = —gipipip, (35a) 



Hic{x) = 9^'^'P 2i{u-d) 
H2{x) = -Sm'tpip, (35c) 



H2c{x) = gSm'tp 



V, (35d) 



H!,,{x) = 5m'^-^^, (35e) 
Hs{x) = -Z^i,^^i,. (35f) 

The Hermitian operator .^J.g^ , which acts on any coordinate space function f(x) standing to the right of it, turns, 
in momentum space, to the factor where k is the momentum conjugated to the coordinate x. The formal 

definition of jf^^Tff^ in coordinate space is given in Appendix A 1. The physical mass of the constituent fermion 

(boson) is m (/tx). The terms H\{x) and Hic{x) in Eq. (34) are the tree level interactions: the usual elementary 
FFb fermion-boson interaction and the contact term on the fermion line between two elementary interaction 
vertices. The terms H2{x), H2c{x), and H2c{x) are associated with the mass counterterm we mentioned above. 

The last, w-dependent, term Hs{x) plays a particular role. Actually, we postulate its form on the basis 
of the observations made in Sec. III. It is needed to eliminate nonphysieal dependence of calculated observables 
on the light front position and may come both from the Fock space truncation and from cutting off divergent 
integrals. A strict derivation of Hz{x) from the "first principles" is however beyond the scope of the present 
article and should be a subject of another work. 

We emphasize that 5ra (but not Sm' !) is the difference between the dressed and the bare masses [see 
Eq. (6) above]. Therefore, it represents the shift of mass (from the bare to the dressed one) due to interaction. 
Since, in gauge theory, the corresponding equations of motion are gauge invariant, the value of 5m defined in 
this way is gauge invariant too, though, generally speaking, the counterterms are not obligatory gauge invariant. 
The gauge invariance of 5m will be analyzed in our calculations. 

We shall apply the Hamiltonian (34) in order to find the composite fermion state vector in the two-body 
approximation. For this purpose, we must know the matrix elements of H"^*{x) between \F) and \Fb) states. 
In practice, it is convenient to represent these matrix elements graphically by means of CLFD diagrams, using 
for calculation of their amplitudes the graph technique rules [3] . According to the latters, matrix elements of an 



^The notation 5m^ means {6m)^ and not 6{m^ 
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interaction Hamiltonian are given by amplitudes of the corresponding diagrams taken with the opposite sign. 
In the following we will refer to these amplitudes as interaction vertices. 

The interaction vertices for the Hamiltonian (34) arc shown graphically in Figs. 3 8 In each diagram the 
four-vector p is the total four-momentum of the initial state. As previously, the solid and wavy lines correspond 
to fermion and boson states, respectively. The dashed lines correspond to spurions — fictitious particles ensuring 
four- momentum conservation at each vertex [3]. There exist three different "elementary" vertices: a three-point 
vertex proportional to g and two two-point ones proportional to 5m and Z^. We will denote these vertices 
by a small full blob, a cross and a cross in an open circle, respectively. The contact interaction between two 
elementary vertices, proportional to — 2(^-k) momentum space, is denoted by a fermion line with a large full 
blob (here k is the four- momentum corresponding to the line going through the large full blob). 

The interaction vertices of -ffi(x), Eq. (35a), are shown in Fig. 3. They describe transitions \F) — > \Fh) 
and \Fb) \F). The term Hic{x), Eq. (35b), in our truncated Fock space, conserves the number of particles. 
It gives rise to transitions \F) \F) and \Fb) \Fb) only. The corresponding interaction vertices are shown 
in Fig. 4. Note that the amplitude of the diagram (c) is infinite since it is proportional to a divergent integral 
over internal boson momentum. The term H2{x), Eq. (35c), produces the only matrix element between one- 
fermion states, as shown in Fig. 5. The term H2c{x), Eq. (35d), has four matrix elements indicated in Fig. 6, 
for transitions \Fb) \F) and \F) ^\Fb). The terms H^J^), Eq. (35e), and Ha^x), Eq. (35f), do not contain 
bosonic field operators and produce only one matrix element each (Figs. 7 and 8, respectively). 

To avoid any misunderstanding, we emphasize that, in agreement with the rules of the graph technique 
for CLFD diagrams [3], the lines with large full blobs representing contact interaction do not correspond to 
propagators, but should be considered as some "complex" vertices, being the products of the factors standing 
at "elementary" vertices. For example, such a line in Fig. 4(c) is given by 



2uj-(jp-k) 



(36) 



where p — k is just the momentum going through the blob. This momentum is completely defined by the 
energy-momcntiim conservation at the vertices, inchiding the spurion four-momentum. 

It is worth mentioning that the spin structures of the matrix elements corresponding to the diagrams 
shown in Figs. 4(c) and 8 are exactly the same. Indeed, these matrix elements (denote them by {H\c) and (.ffs), 
respectively) are given by the following analytical expressions: 

{Hu)=9^Cf,u''^u, {H^) = -Z^u''^u, (37) 
u)-p uj-p 

where 

= 16^ / 2^' 

uj-p 

as follows from Eq. (16), and u (u') is the initial (final) fermion bispinor. The constant Cjc is calculated in an 
explicit form in Appendix B 1. Since Hic{x) and H^(x) come into the full Hamiltonian (34) as a sum, the same 
property takes place for their matrix elements shown in Figs. 4(c) and 8. If we define a new constant 

Z = Z^- g^Cfc, (39) 

the total contribution of (i?ic) and {H3) becomes 

{Hu) + {H3) = -Z u'^^u. (40) 
(j-p 

So, the divergent bosonic loop in Fig. 4(c) is effectively "swallowed" by the counterterm proportional to Z. 



■'Note that two additional contributions should be added in Fig. 4 if higher Fock states are considered. They correspond 
to the case where the both boson lines originate either from the initial or final state. 
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B. Equations in the \F) + \Fb) approximation 



In truncated Fock space we consider in this study the state vector (paip) includes two sectors only. Its 
Fock decomposition has thus the form 



X ^'^'(pi -p- UJTl)2{u-p) 



+ 



(27r)3/2y/2i;: 

(27r)3/^ 



I <^2,aa'(A:i,/c2,P,a;r2)a^,(ki)ct(k2)|0) 



where and ct are the creation operators of the free fermion and boson, respectively, and e^^ = -\/pf + m^, 
£fei = v^k^ + m|, with the appropriate masses. The quantities (\)\ and (/)2 are just the light-front one- and 
two-body Fock components (wave functions), N is the normalization factor. 
The normalization condition for the state vector has the form 

'\>W)'PM = (2^)3(5..' 2£p5(3)(p - p'), (42) 

that gives for the wave functions [3]: 

X^?^la'5(Pl'P''^'^l)'^l,<^s(Pl>P>'^n) + T^^X^ / (t>l^^'ai^l,k2,P,U}T2)^2,cT9{kl,k2,P,UlT2) 

x,5(4)(^i +k2-p- u;T2)2{u-p)dT2^^ = Nd„^,. (43) 

2£fei ZEfej 

The state vector corresponds to a total spin of 1/2 and satisfies the following equation [12]: 

2{u-p) J W-\ujt)^Mp) = - [{Py - M^] Mp), (44) 

where is the free momentum operator, and 

H"*(a;T) = J iJ"*(:r)e-*™'"d4x. (45) 

The interaction Hamiltonian H™*{x) is given by Eq. (34). The mass of the composite fermion is denoted by 
M. In the end of the calculation, we shall take the limit M m. 

Following Ref. [12], we define the one- and two-body vertex functions Fi and r2 by 

Ua'{pi)^iUcr{p) = 2(a;-p)Ti(/)i,^a'(pi,p,a;Ti), (46a) 

Ua'{ki)T2Ua{p) = 2{ui-p)T24>2,acr'{kl,k2,P,U}T2). (46b) 

Due to the delta-functions coming into the decomposition (41), we have pi = p + wti and ki + k2 — p + C0T2, 
that gives 

- M2 (fci + k2f - M2 
= — o 1 '^2 = 7: • (47) 



Now Eq. (44) can be rewritten as 

Q{p) = - j #^"*(u;t) T~^g{p) J, 



(48) 
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where Q{p) is defined by a decomposition analogous to Eq. (41), with the wave functions replaced by the vertex 
functions according to Eqs. (46a), (46b). The operator acting on Q{p) multiplies its one-body part by 1/ri 
and its two-body part by 1/t2. 

The system of equations for the vertex functions Fi and r2 is represented graphically in Fig. 9. The 
loop integrals contain r2 in the integrands. As we shall check in the end of the calculation, r2 is a constant 
for our Fock space truncation, so that we can extract it from the integrals. The latters are thus reduced to the 
on-mass-shell self-energy E given by Eq. (15) at p"^ = M^. Since we should ultimately take the limit M ^ m, 
to calculate E we may put in Eq. (15) p"^ = w? at once. Finally we arrive at the following system of matrix 
equations: 

u{pi)Ti u{p) = u{pi)g{-f,) T2 u{p) 

+ U{p,) Z—^ 2 Ml ^l^P 

+ u{p^)g5m' (-^) (-E)r2«(p) 



u{ki)T2u{p) = u{ki)g ^^t'^^^^ ^^P) 

+ u{k,)g^ (-^) (^^)r2"b) 

(A) ^±^r„(,, po) 

In these equations p^ = M^. For convenience, we have denoted by E the fermion self-energy amputated from 
the coupling constant: 

E = g^t. (51) 

Each term in Eqs. (49) and (50) is in one-to-one correspondence with the graphs in Fig. 9. These graphs include 
various interaction vertices shown in Figs. 3-8, but not all of them. Firstly, we disregarded the graphs with the 
vertex shown in Fig. 4(b), since they contain a three-body (fermion plus two bosons) intermediate state which 
is beyond our two-body approximation. Secondly, by similar reasons, we neglected the graphs with the vertices 
shown in Figs. 6(c,d). Indeed, as already explained in Ref. [12], the counterterm represented by a line with a 
cross is nothing else than a correction to the fermion self-energy indicated in Fig. 2(a). The self-energy involves 
the two-body intermediate state. Therefore, if we have simultaneously one more boson, as shown, for instance, 
in Figs. 6(c,d), such a graph can be considered as a correction to a three-body contribution. Since the latter is 
out of the two-body approximation, we should not take into account the graphs in Figs. 6(c,d). In the system 
of equations (49) and (50) we also combined the contributions from graphs including the interaction vertices 
shown in Figs. 4(c) and 8 into a single term proportional to Zm ijj j (uj-p), according to the discussion in the end 
of Sec. IV A. For simplicity we kept for this term the same graphical notation (a cross in an open circle) as for 
the initial counterterm with Z^. 

Using Eq. (15), the fermion on-mass-shell self-energy E can be represented as 

E= E(p)|^._. =5^A + 5^B^+5^C^, (52) 

TTl UJ'J) 

where, as previously, A = A{w?), B = B^m?), and C = C{m'^). We display explicitly the coupling constant in 
the decomposition of the self-energy. Note that in the general decomposition (52) the term i/j^ is also possible, 
but it does not arise when the self-energy is given by Eq. (15), i.e., when it is linear in the Dirac matrices. The 
coefficients A, B, and C are calculated in Appendix B 1. 
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C. Solution 



The one-body vertex function Fi is proportional to a unit matrix. To avoid the pole terms in Eqs. (49) 
and (50), it is more appropriate to introduce, instead of Fi, the quantity ai defined by 

«i = ^^^2^ (53) 

and related to the one-body wave function by aiUo-' {pi)ucr{p) — (/)i,o-o-' wri). 

Substituting Fi from Eq. (53) into Eqs. (49) and (50) and taking the limit M — > m, we get (the fermion 
spin indices are omitted for shortness) 

= u{pi)g (l - ^Sm'^ (-S) F2 u{p) 
+ u(pi)6m' (l ^(5m' ) 2ma\u{p) 

+ u{pi) ( ) 2maiu{p), (54) 

V i^-P J 



u{ki)V2u{p) = u{ki)g ^1 — ^ ^ 5m!^ 2maiu{p) 



+ ^ik^)9'[-^^i-^)'^Mp)- (55) 

Prom these homogeneous equations we should find 6m' , Z, and (up to a common factor) ai and F2. In contrast 
to Fi, the two-body vertex function F2 has non-trivial spin structure. We represent it in terms of scalar 
components: 

TTl 

u(k\)V2u{p) = biu(ki)u(p) + 62 — u(ki) (/)u(p), (56) 

(j-p 

where 61,2 are some constants to be determined. All other possible spin structures are reduced to those listed 
above. The counting rule which gives the total number N of independent scalar components in the two-body 
vertex is A'^ = (2si -|- l)(2s2 + l)/2 = 2 x 2/2 = 2. The structure u{pi) i/)u{p) for the one-body Fock component 
is not independent since 

U{pi) l/ju{p) = —U{px) l/)TiU{p) = —U{pi){^i- f)u{p) = U{pi)u{p). 

T\ Ti m + M 

We can easily eliminate the spinors in Eqs. (54) and (55) by the replacement u{k) . . . u{l) — > {(i + m) . . .{/ 

I + m), for any momenta k and I. This amounts to multiply both sides by u(k) to the left and by u{l) to 
the right, and then sum over spin projections. After that we take the trace of each equation. In order to get 
the third equation, we multiply the second one by ^m/(w-p) and calculate the trace again. The results are 
expressed through scalar products of the four-vectors entering these equations. In the limit M ^ m, the scalar 
products are given by 

p.pi=m^ p-ki=m'^ - i2'^/2 + {l-x){s-m'^)/2, (57) 

with 

s = {ki+k2f = {p + ujT2f =m'^ + 2{uj-p)T2, a;=— , (58) 

oj-p 

where ki and k2 are the momenta of the two-body system constituents, i. e. the arguments of F2 = 
1^2(^1 ) k2,p,u}T2). Note that because of the different conservations laws the fermion momenta pi (which comes 
into the one-body Fock component) and ki (which comes into the two-body Fock component) have different 
scalar products with p. 
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The system of three homogeneous equations for the three unknown constants ai , bi, and 62 we thus get 
is adduced in Appendix C 1. The determinant of this system is 



Det = -4m^{m^ - sx){l - x) jjm'^ - 2(m + g'^A)5m' 

-g^WB + 2g^C - 2Z) - g^A^] + 2m{g^A + g^B + g^C - Z)] . (59) 

In order the system possesses nontrivial solutions, its determinant must be zero. Equating Det to zero, we 
obtain a quadratic equation for 5m' . It has the following solution: 



5m' = g'^A + m- ^J(m - g^B) \m - g'^B - 2{g^C - Z)]. (60) 

For the other solution (with +yCTT), 5m' does not tend to zero when g ^ 0. We choose therefore the solu- 
tion (60). Substituting this 5m' into the system of equations (C1)-(C3), we can express 61 and 62 through 
ar. 



( I 2(q^C — Z)\ 
bi = 2gmai, 62 = -^mai I 1 - W 1 ^ _ j • (^^) 

We see that both bi and 62 do not depend on s and x. This justifies the extraction of r2 from the integrals 

when we derived Eqs. (54) and (55). The solutions (60) and (61) arc nonpcrturbativc, since their expansions 
in powers of the coupling constant g contain an infinite number of terms. Eqs. (54) and (55) are valid for any 
g (for the given Fock space truncation) provided the expressions under the square roots are non-negative. 

According to Eqs. (46b) and (56), the components which determine the two-body wave fimction are 
obtained by dividing bi and 62 by the factor 2{uj-p)t2 = (s — m^). As a result, the wave function depends on w 
through the term 



62 m 



(fci) i/juip). 



s — m? uj-p 

This dependence must disappear in any physical observable, in particular, in the residue of the wave function 
at the pole s = m^. Since, as follows from the second of Eqs. (61), 62 is a constant, the only way to ensure this 
property is to require 62 = 0. From here we find the constant Z = g^C and 

5m'=g^{A + B). (62) 

Taking into account Eqs. (33) and (39) yields the final solution 

61 = 2gmai, 62 = 0, (63) 



5m = g\A + B + C + Cfc), Z^=g\C + Cfc). (64) 

The three constants ai, 5i_2 determine the one- and two-body vertex functions Fi and The components 
01, (j)2 of the state vector (41) are related to Fi, F2 by Eqs. (46a) and (46b). Explicit expressions for 5m and 
Zi^ are given in Appendix B 1 . 

Substituting the solution for the wave functions into the normalization condition (43), we find the nor- 
malization factor N: 

.111. 9^ f Tr[iUi+m)ii + m)\ ^t,-.,, , ,, d^kid^k^'] 

^ = ^™ ^H'^^y {s-m'f h ^Hk.+k2-p-.r2)2i..p)dr2^^^^^ . (65) 

The integral in Eq. (65) diverges logarithmically. Explicit formulas for N are given in Appendix B 1. Note that 
since the factor a/JV is proportional to ai, as well as each of the wave functions, the state vector (41) does not 
contain the constant ai at all. So, the final solution for the state vector is completely determined by the values 
of the coupling constant and particle masses. 



14 



We remind that we did not renormalize the fermion field by introducing the factor Z2 as in the first of 
Eqs. (2). The fermion field renormalization is not needed in our approach, because the factor N does the same 
job for the whole state vector. The two factors N and Z2 are tightly related to each other. Namely, let us define 
the constant Z2 as 



Z2 = 



N 



It turns out that it is expressed through the self-energy in the standard way: 

d 



Indeed, using the representation (15) of the self-energy, we get 



2m 



+ 2m 



m 



(66) 



(67) 



(68) 



We used the fact (proved in Appendix B 1) that C(p^) does not depend on and that = ^j^p = 0. 

Substituting the expressions (B2a) and (B2b) for »4(p^) and H(p^) into Eq. (68), calculating the corresponding 
integrals and comparing the result with Eq. (B7) for N, we find the relation 



N 



= 1 



_d_ 



S(P) 



(69) 



From Eqs. (66) and (69) immediately follows Eq. (67). Note that the latter is just a standard expression for 
the traditional counterterm (see, e. g., Ref. [18]). 

In perturbation theory, taking into account that E(p) is of order ^ we can transform Eq. (67) to 



^. = l+^E(p) 



1 r , ,aM(2)(p) 



(70) 



where A^^(p) is defined by Eq. (14). The last equality follows from the fact that OS jc{p)/d ^ = 0. The r. h. s. 
of Eq. (70) coincides with that of Eq. (30) obtained previously in perturbative renormalization technique. 

The solution (60), before specifying the value of Z, is true nonperturbative, whereas the final solutions (62) 
and (64) coincide with the perturbative ones. It is not a surprise since the self-consistent calculation in the 
\F) + \Fh) approximation generates all contributions to 5m. (or 5m') of order g^, and nothing else. A real 
difference would appear when higher Pock sectors are taken into account. 

The value (62) of 8m' is in agreement with the solution given in Ref. [14] . 



V. GAUGE BOSON IN THE FEYNMAN GAUGE 



A. Contact terms 



We shall follow here very closely the procedure detailed in the previous section. The light-front interaction 
Hamiltonian for a system of spin 1/2 fermion coupled to gauge bosons is discussed in Appendix A 2. It has the 
general form 



2i{uj-d) 



—Sm'ipil) + gSm'ip 

-Z^%p ip 
+H'{x). 



^ ip 



2i{oj-d) 2i{uj-d) 



2i{u)-o) 



(71) 
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and differs from the Hamiltonian (34) by the replacement ^ /p and by the additional contribution H'{x) 
which represents contact terms associated to the gauge field. An explicit form of H'{x) depends on the gauge. 
Generally speaking, H'{x) taken in an arbitrary gauge and expanded in powers of g gives rise to an infinite 
number of contact terms, even in truncated Fock space. 
For the Feynman gauge, H'{x) acquires the form 



H'p{x) = H[p{x) + H'^p{x) + H'^p{x) 



(72) 



with 



2i{uj-D) 



H2f{x) = gSm'ljj 







2i{LU-D) 



V', 



2i{uj-D) 2i{uj-D) 



V', 



(73a) 
(73b) 
(73c) 



where the Hermitian operator 



i(ui-D) 



is defined in Appendix A 2 a. An expansion of Hp {x) in powers of g starts 



with a term of order g (we took into account that the corresponding expansions of 5m' and of the operator 
begin with terms of order g^ and g, respectively). In full (non-truncated) Fock space this series does not 
break and the number of contact terms is infinite. However, it becomes finite in truncated Fock space. Since we 
are working within one- and two-body sectors, we have to know the non-zero matrix elements of H'p{x) between 
these sectors. Let us consider first the term H[p{x). We can represent H[p{x) as the following series in powers 
oig: 



H[p{x) 



711=0,712=0,711+712)0 



E 



,"l+"2 



Hn 



(74) 



tjj ijp 



UJ-Lp 



iuj-d 



iuj-d 



LO-Lp 



(75) 



Each term Hn^n^ contains ni + na + 2 boson field operators ^pi,. When we calculate matrix elements between 
states with the fixed number of particles, these operators must be contracted with external boson operators or 
among themselves. It is easy to see that contractions between two different operators lo-p or between lo-p and 
ifp give zero in the Feynman gauge. Indeed, since (0|[</5ju, (pi,]|0) ~ —g^iv, the "internal" contractions between 
the operators w-tp in Eq. (75) are proportional to cj^^iv^g^^, that is zero. Contractions between w-ip and jp also 
disappear because they are proportional to 7^ ^cj'^g^jy =^ ^ = = 0. The only remaining possibility is to 
contract the operators Lu-p with the external boson field operators. The number of the latters is fixed and docs 
not exceed two in our approximation. Hence, we must have ni + n2 < 2. This just means a truncation of the 
series (74) to a polynomial of the fourth order in g. Evidently, the same arguments hold true for each of the 
terms H2p{x) and H.^p{x), and, hence, for the operator H'p{x) as a whole. Note that H'p{x) does not have 
non-zero matrix elements within the one-body sector, because it contains at least one operator Lu-ip. 

Considering matrix elements of the operator H'p(x) between the states \F) and \Fh), we find that each 
piece H[p{x), H'2p{x), or H^p{x) produces eight different interaction vertices. For example, in Fig. 10 all 
possible interaction vertices generated by the piece H[p{x) are shown. Each "elementary" vertex (small black 
blob) corresponds to the factor gj'^, the propagator of boson with four-momentum q being —g^jy6{q^)9{tjj-q). The 
other notations are the same as those adopted for the scalar boson case, except for new factors indicated by large 
open circles with two fermion and one boson legs. These factors proceed from different powers of the operators 
in Eq. (75). A large open circle corresponds to the factor ±g^^, where k is the four-momentum of the 
boson line connected with this open circle. The sign plus (minus) is taken for an outgoing (incoming) boson 
line. Additionally, each matrix element must be multiplied by a combinatorial factor l/(ni!n2!), where ni (722) 
is the number of large open circles standing to the left (to the right) from the large full blob. We have omitted 
two diagrams which differ from those shown in Figs. 10(g) and (h) by the interchange of the initial and final 
boson lines because such an interchange does not affect analytical expressions for the corresponding amplitudes. 
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So, to incorporate the contributions from the omitted diagrams it is enough to multiply the amplitudes of the 
diagrams in Figs. 10(g), (h) by a factor of two, just cancelling the combinatorial factor 1/2! 

As far as the interaction vertices produced by the terms H2p{x) and H'^p(x) are concerned, they can be 
easily obtained from the diagrams of Fig. 10. Indeed, to get the vertices for H2p{x) one should substitute in 
all these diagrams one of small full blobs by a cross and transform the boson loop to a free boson line incoming 
to (or outgoing from) the remaining small full blob. In such a way each of the diagrams in Figs. 10(a) (d) 
transforms into two new ones, while the diagrams shown in Figs. 10(e)-(h) give nothing because all of them 
contain after transformation three boson legs. The interaction vertices for H^p{x) are obtained from those in 
Fig. 10 simply by substituting all small full blobs by crosses and by removing the boson loops. 

Note that similarly to the scalar boson case, any line with a large (full or open) blob, or with two of 
them, is not a propagator, but only a product of the factors corresponding to the elementary vertices (large and 
small full blobs, and large open circles). For example, the internal solid line in Fig. 10(a) corresponds to the 
analytical expression 

where p = pi + ki — loti = p[ — u!t[ is the total momentum of the initial or final state. To calculate the 
full matrix element, one should sandwich Eq. (76) with bispinors. then multiply it by the boson propagator 
—9niyS{q^)0{co-q) and by the polarization vector e^(A;i), integrate over d'^q/{2-K)^ and multiply the result by 
(-!)■ 

It is important to note that although we have calculated matrix elements of H'p{x) between the states 
\F) and \Fb), any of the interaction vertices obtained from H'p{x) includes intermediate states with at least two 
bosons or with one boson and a cross. Repeating the same arguments as those exposed for the scalar boson 
case, we have thus to conclude that the whole contribution from H'p{x) is beyond our two-body approximation. 







2uj-{p - q) 



B. Equations in the \F) + \Fb) approximation 

The Fock decomposition of the state vector is similar to the one for scalar bosons, given by Eq. (41). 
The only difference is the dependence of the boson creation operator and the Fock components on the boson 
spin index A. As usual, we define the one- and two-body vertex functions Fi and Fj by [cf. with Eqs. (46a) 
and (46b)] 

Ua'{pi)Tiu„{p) = 2{uj-p)Ti(pi^^^i{pi,p,ujTi), (77a) 
Ua'{ki)T^e^^{k2)ua{p) = 2{co-p)T2(pl^^,ikuk2,p,ujT2). (77b) 

We should now establish, in the two-body approximation, the whole set of matrix elements of the Hamil- 

tonian (71) between the one- and two-body sectors. Since the Hamiltonian (71), excepting for the term Hp{x), 
looks quite similar to that for the scalar boson case, Eq. (34) , we can make use of this analogy. So, we can take 
directly the graphs shown in Figs. 3-5, 6(a,b), 7, 8, substituting g — > 57'' everywhere, and add to this set the 
contributions from Hp{x), indicated in Fig. 10, which are compatible with our Fock space truncation. However, 
as was noticed in the end of Sec. V A, these contributions correspond in fact to many-body sectors and should 
therefore be disgarded. 

For further calculations we will need the fermion self-energy. In the Feynman gauge it is 

{-Qi^i^h^ifl- ^+ + m)Y e[uj-{p - k)] d^k 



(27r)3 J 2w(p - k)T 2ek 

= 1' K.) r = 9' A + + g^C^ (78) 

with r defined by Eq. (17). Note that the constants A, B, and C here do not coincide with those for the scalar 
boson case, defined by Eq. (52). From Eq. (78) follows 

= -9,.i9'A, + g'B.l + 9'C,^) ^ g^%^^ (79) 
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where 



A = -AAx, B = 2Bi, C = 2Cx. (80) 

The constants A, B, and C are given by Eqs. (B12a)-(B12c) in Appendix B2. 

We can now write down a system of equations for the vertex functions. The graphical representation 
of these equations is shown in Fig. 11. It is very similar to the one for scalar bosons, as shown in Fig. 9. 

The only formal difference consists in the appearance of fonr-dimcnsional indices in boson propagators and 
vertices. Generally speaking, the vertices Fi and Ti^ do have momentum dependence. Because of our Fock 
space truncation, they are however constants. At the present stage we take this property as an ansatz, and it 
will be justified by the final results. 

In an analytical form, this system of equations reads 

u{pi)T,u(p) = M(pi)57'^(-E^^)r^u(p) 



+ ^^^^^ V ^ ) m^-A'P ^^^(^^ 



+ u{p,)Sm [-^j—^-^T,u{p), (81) 



(^+ i^Ti + m) 



+ u{k,)g^j^e^^{k2) [-^) Y{-K)^Hp) 

+ Uik.)gSmYe^,ik.) (-^) ^^^^^r.«(,). (82) 

As explained in Sec. IV A, we absorbed the (divergent) term proportional to Z^^ and the contribution from 
the diagram shown in Fig. 4(c) (with the scalar boson replaced by the gauge one) in the unique counterterm 
proportional to Z, defined by Eq. (39), where C/c is now given by 



^ cj-p ~ (2^)3 J 2Loip - k) 2ek ■ ^^^^ 
The constant C/c is calculated in Appendix B2 [see Eq. (B12d)]. 

C. Solution 

In order to solve Eqs. (81) and (82), it is convenient to decompose the vertex functions in invariant 
amplitudes, similarly to what was done in Sec. IV. The structure of u{pi)Tiu{p) coincides with the one for the 
spinor-scalar interaction [cf. with Eq. (53)]: 

-^^^^^=a,u{p,)u{p), (84) 

where a\ is a constant and the factor [m? — M^)~^ reflects the behaviour of Fi in the limit M ^ m. The 
number of independent invariant amplitudes for the vertex function u{ki)\r2e^{k2)]u{p) coincides with that 
for the reaction 1/2 + 1/2 + 1. However, one should take into account that in the Feynman gauge the 
vector boson wave function has Jour independent components. So, the total number of invariant amplitudes is 
(2 X 2 X 4)/2 = 8. We choose the following set of invariant amplitudes: 
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617'' + 62 \-b; 

ui-p 



u>-p 



+ b. 



+ b5— +be— +b7- 

m m co-p 



+ bs 



w-p 



e^(fe)w(p). 



(85) 



One can easily check that in our approximation the structures proportional to p^ and fcf do not contribute at 
all. Indeed, whatever term from Eq. (85) we substitute into Eq. (82), the result will always include structures 
proportional to 7'^, w'^, (Z'7'', or j^w'' only, i. e. the four structures proportional to 65-8 must be absent. Although 
all the four coefficients 65-8 can be obtained automatically, as a solution of a system of four linear equations, 
with the result 65 = feg = ^7 = ^8 = 0, we may set b^-% = from the very beginning in order to simplify the 
calculations, and reduce Eq. (85) to 



tZ(fci)[r^e^(fc2)]«(p)=«(fci) 



617'' + b2 



uj-p 



+ &3 +64 



ijj-p 



iv-pf 



el{k2)u{p). 



(86) 



To solve the equations (81) and (82), we repeat the same steps as done in Sec. IV C for scalar bosons. We insert 
Eqs. (84), (86) into Eqs. (81), (82), take the limit M — > m, and substitute the spinors and polarization vectors 
according to u{pi) . ..u{p) {pi + m) . . .{p + m), u{ki) . ..u{p) (/;i + m) . . . (j6 + to), ^^^(fca) -.9^^'- 
After these transformations we take the trace from Eq. (81). Then we multiply the transformed Eq. (82) in 
turn by 7^^ , muj^ /("^'P), ™7^ f^li^'P)-, and p^ /m, and take traces again. As a result, we get a system of 
five homogeneous linear equations for the five coeflicicnts ai, 61-4, which is given in Appendix C2. For finding 
nontrivial solutions of this system, we have to demand its determinant to be zero: 



Det = Am^irri^ - sx)x{l - x^) |(5to'^ - 2(m - Ag^Ax)5m' 

-4g^[Biig^Bi + 2g^Ci - Z) - Ag^Af] - 2m{Ag''Ai - 2g^Bi - 2g^Ci + Z)} = 0. (87) 

Solving Eq. (87) with respect to Sm' , we find two solutions. Evidently, we should choose the solution which 
gives 5m' — > as g — > 0. Expressing Ai, Bi, and Ci through the constants A, B, and C by means of Eqs. (80), 
we get 

Sm' = -Ag^Ai +m- ^/{m - 2g'^Bi)[m - 2g'^Bi - 2{2g'^Ci - Z)] 

= g'^A + m- V(m - g'^B) [m - g'^B - 2{g'^C - Z)]. (88) 

Substituting Eq. (88) into Eqs. (C4)-(C8), we obtain 

61 = 2gTOai, (89) 



62 = -263 = -2gmai [l-^l- ^^^^-^j , 64 = 0. (90) 

Like to the scalar boson case, we require the residue of the two-body wave function at the pole s = to^ to be 
w-independent. This leads to the conditions 62 = 63 = 64 = which can be satisfied if Z = g^C. We thus arrive 
to the following final solution for the vertex functions and the two counterterms: 

61 = 2gmai, 62-8 = 0, (91) 



6m = g^A + B + C + Cfe), Z^=g\C + Cf,). (92) 

Eqs. (92) formally coincide with Eqs. (64) foimd for the scalar boson case. However the constants A, B, C, and 
Cfc theirselves depend on the boson type. Due to this reason the exact values of the counterterms 6m and Z^^ 
here are different from those in Eqs. (64). They are given in Appendix B 2. 



19 



VI. GAUGE BOSON IN THE LIGHT-CONE GAUGE 



A. Contact terms 



The interaction Hamiltonian has the same form (71) as in the Feynman gauge, except for the gauge 
dependent term H'{x). In the LC gauge where w-</p = we have (see Appendix A 2 b) 

H'lc{^) = g\i^ (^V') 2(-^(Vi ^^)- (93) 

In contrast to all other gauges, this contribution to the light-front Hamiltonian generates only one contact term, 
of order g^, associated to the propagator of the gauge field. 

As far as the non-zero matrix elements of H'j^q{x) are concerned, the situation here is simpler than in 
the case of the Feynman gauge. This part of the Hamiltonian conserves the number of particles in truncated 
Fock space which we consider here. It generates two interaction vertices in the one-body sector, as shown 
diagrammatically in Fig. 12. 

The fcrmion contact term with the boson loop, shown in Fig. 4(c) (with the scalar boson replaced by the 
gauge one) is now 

where 



^ iv-p ~ {2n f J 2a;-(p-fc) 2£fc ' ^^^^ 



d^,u = -g^v + — ; (95) 

results from the vector boson propagator in the LC gauge. 

The sum of the two (divergent) fermionic loops shown in Fig. 12 represents a contact correction, on the 
boson line, to the self-energy contribution. It can be written as 

2„ m i/j _ g"^ f f i/} (( i/j i/j (/i \ (Pk 



^ iv-p - (27r)3 J [2[u;-{p - k)]^ 2[uip + k)]^ 2sk ' ^^^^ 

The expressions (94) and (96) have the same matrix structure as the one coming from the cj-dependent coun- 
terterm proportional to Z^. We can therefore again absorb them, as in Eq. (39), into a unique term which, for 
simplicity, we also call Z: 

Z = Z^- g^Cj, - g^Cbo. (97) 



B. Equations in the \F) + \Fb) approximation 

The contribution (96) should be excluded from the equation for the last Fock component The reason 
for this is that H'^(~,{x) defined by Eq. (93) describes corrections to bosonic propagators. These corrections to 
r2 can be regarded as " foot- prints" of contributions from the three-body sector containing one fermion and two 
bosons, which is outside of our model space. 

The fermion self-energy has now the form 

(27r)3 J 2w{j> - k)T 2£k 

^ r^^^^Y = g'A + g'sl + g'c"^, (98) 

where d^i, and r are given by Eqs. (95) and (17), respectively. 
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The symmetric tensor d^j^ is transversal to both uj and the boson momentum fc. As a consequence, the 
operator Y,^^ in Eq. (98) is a symmetric tensor transversal to uj. Besides that, it contains only zeroth and first 
powers of 7-matrices. Due to all these, the tensor has the following general structure: 



= - ( 9tiv - , _ ) ( 9 M + g^Bi — ] + 9 Ci— {g^^ ^ - 7^w^ - 7^0;^) 



m 



+5 C'l — 1/)^ - 7^0;^ - 7^0;^ + — ^ g^A2 + g^B2— + g^C2 , (99) 

determined by seven scalar functions Ai, Bi, Ci, €[, A2, B2, and C2. As usual, we define — g^Y/^ . 
Multiplying Eq. (99) by 7^ to the left and by to the right, contracting over the indices /x and and 
comparing the result with the r. h. s. of Eq. (98), we get 

A = -2Ai, B = 2Bi, C ^2{Bi+B2-5C\-'iC[). (100) 

The values of A, B, and C are given in Appendix B 3 by Eqs. (B16a)-(B16c). 

We thus arrive, in the LC gauge, at a system of equations for the vertex functions, which is quite similar 
to that for the case of gauge bosons in the Feynman gauge, shown graphically in Fig. 11. Because H'^q{x) does 
not produce any matrix elements of new structure as compared with the case of the Feynman gauge, we can 
get the analytical form of the equations for the functions Fi and r2 directly from Eqs. (81) and (82). In the 
LC gauge the system of corresponding equations reads 

«(pi)Fi«(p) = u{p^)gY{-^i^)n<p) 

+ u{pi)5m — — Fiu p 

— 

I -( \ ( ry'm. iZiri + to) 

+ m(Pi) z ^ 772-^riu(p) 

\ uj-p J — M-' 

+ u{p,)gSm' (-^) n-Y^^nuip) 

-^^^^^A-d^Y-^^^ ''''' 

u{h)[T^,e'^{k2)]u{p) = n(fci)fl7^e^(fc2) ^^^t' j^^^ ^Mp) 

+ u{k,)g^rel{k2) (-^) 7''(-S^,^)r^i(p) 

+ .(.O.W7'^e^(^2) (-^) ^^it^^r^ ^Mp)- (102) 
Note however that the polarization vectors e^{k2) here differ from those in the system (81), (82). 

C. Solution 

As above, the structure of the one-body vertex function Fi is given by Eq. (84). The number of inde- 
pendent invariant amplitudes for the vertex function u(fci)[r2e^(A:2)]u(p) can be calculated analogously to the 
case of the Feynman gauge, taking into account, however, that in the LC gauge the vector boson polarization 
vector is orthogonal to both u) and the boson momentum. The conditions e^{k2)-to = and e^{k2)-k2 = leave 
two independent components. The total number of invariant amplitudes is thus (2 x 2 x 2)/2 = 4. We choose 
the following set of invariant amplitudes: 

S(fci)[F^e;i(fc2)]u(p) = u{k^) Ul" + 62^^^^ + ?>3- + 64^^^^) el{k2)u(p). (103) 
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Again, in our approximation, the structures proportional to do not arise at all. We can therefore put 
63 = 64 = and reduce Eq. (103) to 

M(fci)[r^e^(fc2)]«(p) = u{ki) (^617" + b2^^^^^ e^^{k2)u{p). (104) 

We insert Eqs. (84), (104) into Eqs. (101), (102) and replace, as before, the spinors u{pi), u{ki), u{p) by (^i+m), 
{^1 + m), + to), respectively, and the polarization vectors e^{k2) by -[g,j.,i' - [katiOJ^i' + fc2^'t^^)/(t^-fc2)]- 
After that we take the limit M ^ m and calculate the trace from Eq. (101). Then we multiply Eq. (102) by 7** 
and by to ((67^ / {uo-p) successively and take traces again. In calculating the traces we need to know the scalar 
products p ■ k2, ki ■ k2, and lu ■ k2 in addition to those given by Eq. (57). These scalar products can be easily 
expressed through the variables s and x [see Eq. (58)] by using the conservation law kx+k2 = p + uiT and the 
conditions kf = p^ = rn?, A;f = 0, = 0. The result is 

p-k2 = , kvk2 = — - — , = l-x. (105) 

2 2 LO-p 

Finally, we get a system of three homogeneous linear equations for the coefficients ai, 61,2- It is given in 
Appendix C3. 

Equating the determinant of this system to zero, we come to a quadratic equation for Sm': 

Det = ^sm" + 2(2,^ A, - m)Sm' 

+4g^[g^Al - Bi(35'Bi + 2g^B2 - Wg^d - 8g^C[ - Z)] 

-2m{2g'^Ai - Ag'^Bi - 2g^B2 + lOg'^Ci + 8g^C[ + Z) | = (106) 

with the solution (we choose the minus sign at the square root, as explained above) 

Sm' = -2g'^Ai + to - J {m - 2g'^Bi)[m - 2{3g^Bi + 2g^B2 - lOg'^Ci - Sg'^C'-^ - Z)] 



= g^A + m- ^J(m - g'^B) [to - g'^B - 2{g^C - Z)], (107) 

where we used the relations (100) between the coefficients. The substitution of Eq. (107) into Eqs. (C9)-(C11) 
allows to express the functions 61,2 through ai. In our two-body approximation the functions 61,2 are constants: 

bi = 2gmai, (108) 



/ 2(2ff2Ei + 2g^B2 - lOg^C, - 8g^C'^ - Z) 
b2 = gmai\l-Al ^ 



TO — 2g'^Bi 



If we impose the condition Z — g^C = on the renormalization constant Z which is still free, the function 62 
turns to zero, as it should. Taking into account Eqs. (33) and (97), we find the following final solution: 

61 = 2gmai, 62-4 = 0, (110) 

5m = g^{A + B + C + Cf, + Chr), Z^ = g'^{C + Cf, + d^). (Ill) 

It formally differs from the corresponding solution (91) (92) for the case of the Feynman gauge by the presence 
of the boson contact term Cbc- Note however that the values of A, B, C\ and Cjc are different for the two 
gauges. Explicit expressions for Sm and Z^ are given in Appendix B3. 

Though the solution for the state vector in the LC gauge, Eqs. (110), formally coincides with that obtained 
for the Feynman gauge and given by Eqs. (91), the normalization factors of the corresponding state vectors are 
different. These normalization factors are calculated fully analogously to how it was done for the scalar boson 
case. 
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VII. ANALYTICAL RESULTS AND DISCUSSION 



Let us denote the contact term contribution by Cc- 

(J _ i C'/c for scalar boson and for gauge boson in the Feynman gauge, (112) 

~ \ Cfc + Cbc for gauge boson in the LC gauge. ^ ' 

Then the expressions for and 6m in terms of A, B, C, Cc obtain the same form for all the cases discussed 
above: 

5m = g^{A + B + C + Cc), (113) 
Z^^g^{C + C,). (114) 

The constants A, B, C, and Cc are however different for different bosons. For the scalar boson case all the 
constants are given in Appendix B 1. For the cases of gauge boson in the Feynman and LC gauges they are 
given in Appendices B2 and B3; respectively. Substituting their explicit values into Eqs. (113) and (114), we 
find the analytical formulas for Z^^ and 5m. 



For scalar boson: 



„ „^ Q^m 

fa.= = -|^log-j, (115) 



For gauge boson in the Feynman gauge: 



For gauge boson in the LC gauge: 



<5m-^ = $flog^, (119) 



87r2 
32^ 



Z^^ = ^log^. (120) 



Let us now discuss what happens with fermion-boson scattering amplitudes calculated in CLFD within 
the two-body approximation. Although this question, strictly speaking, goes beyond the scope of the present 
paper, we think that our results are rather instructive and may be useful for future investigations. In Sec. Ill 
we considered the perturbative expansion (up to terms of order g^) of such an amplitude for the scalar boson 
case and showed that, after introducing the counterterm Z^^, it became w-independent. Below we perform a 
similar analysis of the scattering amplitudes without making expansions in powers of g. 

In the two-body approximation the fermion-boson scattering amplitude M can be obtained by sandwich- 
ing the fermion 2PGF, V^p), with the corresponding bispinors. The 2PGF, in its turn, is easily derived from 
Eq. (12) by putting there Z2 = I and by adding the counterterm (27) to dm. As a result, we find 

M = igMp'iMPHPi) = -qMp'i) f + ^ ^(p^^^ (121) 



— + 



5m-7U(2)(p)+Z„(^-l)] (^ + m) 



where F = F' = 1 for scalar bosons, F =/(fci), F' =j^.{k[) for gauge bosons, ki and k[ are the boson four- 
momenta before and after scattering. The superscript of the mass operator indicates that we are working in 
the two-body approximation. Using Eq. (112), we can write M^'^\p) in the form common for all boson types 
[cf. with Eq. (18)]: 



M^'Hp) = 9' {a{p') + B{p^)l + [C{p') + CJ^ 



m oj-p 



(122) 
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The coefficients A, B, and C in Eqs. (113) and (114) are, respectively, the values of the functions A{p^), B{p'^), 
and C{p^) at = rr? . In spite of that Eqs. (121) and (122) have the same form for all boson types, these 
functions, as well as the constants (5m, Cc, and Z^^, arc different for different bosons. 

Substituting Eq. (122) into Eq. (121) and using Eqs. (113) and (114), we see that the w-dependence of 
the amplitude is determined by the term 

{g" [C(p2) + Ce] -Z.) !^(^ + m) [C(/)-C(m2)] ^(|( + m) (123) 

in the denominator of Eq. (121). In Appendix B 1 we proved that C(p^) — C(to^) = for the scalar boson case. 
For the case of gauge boson in the Feynman gauge the situation is completely the same, since the function 
C(p^) differs from that for scalar bosons, given by Eq. (B2c), by a factor of 2 only (see the argumentation in 
Appendix B2). So, the quantity (123) equals zero for both scalar bosons and gauge bosons in the Feynman 
gauge, and the scattering amphtude (121) does not depend on uj at all. 

The latter statement is not true for the case of the LC gauge. The function C(p^) can be obtained by 
substituting the integral from Eq. (98), without putting = rn^, into the last of Eqs. (Bl). Technical details 
are completely equivalent to those exposed in Appendix B 1. Thus, we get for the LC gauge: 

l_ kl{x--2x + 2) + m-x--p-.Hl-x? 

^^^>~ Wn^mJo "-Js xHl~x)[kl+m^x-p^x{l-x)] ' ^ ' 

Subtracting from Eq. (124) the value C(m^), we arrive at the expression 

C(v') - C(m') - ~ "^'^ dk^ C dx kl{2x---ix + 2)+m^x^ 

L[p ) L[m)- '^^^ J,'' x[kl+m^x^][kl+m^x-p^x{l-x)y ^ ^' 

The latter integral diverges logarithmically at x = and can not be zero. From here follows that the scattering 
amplitude (121) for gauge bosons in the LC gauge, generally speaking, depends on uj. Such an unpleasant 
feature is however not specific to LFD, but also takes place in 4DF formalism, as explained in Ref. [19]. This 
is caused by the bad regularization of the singularity in w-k to which the function C(p^) is sensitive. 
Now we can summarize the results obtained above. 

• LFD requires an additional counterterm in the Hamiltonian, depending on the light front orientation 
(i. e., on the four-vector w). In the two-body approximation this counterterm in momentum space has the 
form —Z^^. In the nonperturbative calculations presented above it eliminates the w-dependence of the 
two-body wave function residue at the pole s = rr?. 

• We have found, for the truncation \F) -\- the true nonperturbative solutions for the mass counterterm 
(5m and the composite fermion state vector for the cases of scalar and gauge bosons. Expanded in a series 
in powers of the coupling constant g, these solutions contain all degrees of the latter. Besides that, the 
two-body wave functions have a non-zero spin components explicitly depending on lo. The additional 
counterterm, Z^^^ eliminates these components. After that the solutions (115)-(120) coincide with the 
perturbative ones, as they should if the \F) -I- \Fh) Tamm-Dancoff truncation is used. A real difference 
between perturbative and nonperturbative results must appear when higher Fock sectors are taken into 
account. In the latter case one can hardly expect the vertex functions to be constants. Consequently, 
there is no any reason to demand that all their w-dependent components equal zero identically in order 
to get an w-independent residue of the two-body wave function at its pole. 

• The mass shift 5m is gauge invariant, as is seen from the coincidence of Eqs. (117) and (119). 

• The mass shift (5m does not coincide with the corresponding value found in Ref. [6], since it depends on 
the choice of cutoffs. Our cutoffs are different from those used in Ref. [6]. We show in Appendix D that 
our formulas with the cutoffs taken from Ref. [6] reproduce 5m obtained there. 

• The value of the counterterm Z^^ depends on gauge, as expected. 

• In spite of the fact that intermediate calculations involve strongly divergent integrals [see, e. g., Eqs. (B4c), 
(B4d), (B12c), (B12d), (B16c)-(B16e)], the final results for 5m and Z^ contain logarithmic divergences 
only. 
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VIII. CONCLUSION 



We considered a spin 1/2 fermion state consisting of the bare fermion F coupled to bosons b with spin 
either zero or one, in the approximation where the first two Fock sectors (|F) + \Fb)) are kept in the state 
vector. The spin one bosons arc the gauge ones. 

Although the state vector of the composite fermion is simplified, the study of this model is of utmost 
importance in order to settle the various pieces of the light-front Hamiltonian we start from. While the light- 
front Hamiltonian is simple at the tree level, this is not true anymore when higher order processes are considered. 
Due to the explicit covariance of our approach, we exhibit in the present study the exact operator structure 
of the counterterms needed for the renormalization. In the two-body approximation two counterterms are 
required to renormalize the theory: a mass counterterm and a specific one for LFD, which explicitly depends 
on the orientation of the light front plane (i. e. on the four- vector lu). We calculated them in an explicit form. 
While the mass counterterm is analogous to that which appears in standard 4DF approach, the specific LFD 
counterterm has no such an analogue. 

These two coimtcrtcrms generate additional contributions to the light-front Hamiltonian, due to their 
coupling with contact interactions. The latters also depend on the light front plane orientation. The solutions 
we found for the couplings of scalar and gauge bosons to a spin 1/2 fermion are in complete agreement with 
the CLFD perturbation theory, although our calculation does not rely on perturbative expansions. The condi- 
tion imposed on the mentioned above specific counterterm eliminates the w-dependence of the two-body wave 
function residue at its pole. 

We explicitly showed that the mass counterterm 6m is gauge independent. We calculated the state vector 
in two gauges — the Feynman and the LC ones. Though finally the number of the two-body wave function 
spin components turned out to be the same in both gauges, this is a peculiarity of the \F) + \Fb) Fock space 
truncation. Such a property should disappear when higher Fock components are taken into account. At the 
same time, physical observables (e. g, the electromagnetic form factors) must be the same in any gauge. 

The results reported here are encouraging in the perspective of tackling physical composite systems in 
QED and QCD. This will be a subject of forthcoming researches. 
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For non-zero spin fields, the light-front Hamiltonian should incorporate the so-called contact interactions. 
Their origin and derivation are explained, e.g., in Ref. [2] (in the LC gauge, for the gluon field). For completeness, 
we briefly explain this derivation, in the Feynman gauge too. 

We start with the case of scalar boson field. The corresponding field theory Lagrangian, at the tree level, 

is 



where 4* and $ are, respectively, the Heisenberg spinor and scalar field operators. In a standard way we obtain 
the energy-momentum tensor 



ACKNOWLEDGEMENTS 



APPENDIX A: LIGHT-FRONT HAMILTONIANS 



1. Scalar boson 




(Al) 



^ dC dC 



8"^ + 



dC 



(A2) 
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and the four-momentum operator 



r^ = J da^{x)e''P. (A3) 



Integration in Eq. (A3) is performed on the three-dimensional space element orthogonal to the " time" direction. 
The role of time is played in CLFD by the invariant combination uj-x. To simplify the subsequent algebraical 
manipulations it is convenient to introduce temporarily a particular value^ of the four- vector oj, namely, uj = 
(1,0,0,-1), corresponding to the usual form of noncovariant LFD on the plane x° + x^ = 0. In the final 
results we will return to explicitly covariant notations. This is simply achieved by the replacement of the 
plus-component of any four-vector a by the contraction Lu-a. 
Let us introduce also the light-front coordinates 

x+ = x°+x^, x-=x°-x^, yi^ = {x\x'^}. (A4) 

The scalar product is defined as 

a^b" = ^a+b- + ^a-b+ -a^ -h^. (A5) 

The metric tensor is 

g+-=g-+ = 2 g+-=g-+ = ^ g^^ = g-'^ = g^^ = g22 = -1, (A6) 

all the other components being zeroes. In such a formalism the and 'P"'" components of the momentum 
operator are free, while the V~ component contains interaction. "P" is given by the expression 

V- = ^J <fx^dx-Q+-. (A7) 

Prom Eqs. (A2) and (A3) with the Lagrangian (Al) we get 

V~ = Sx^dx- -i*7-a+* -f- 2^*7"^ • a^* + 2m** -|- {d^^f + /x^^^ - 2g^^^^ . (A8) 

The time evolution of the operators * and $ is governed by the equations of motion obtained from the 
Lagrangian (Al). Thus we have the following Dirac-type equation for *: 

(i ^ - m)* = (A9) 

The spinor field * has four components. However, it turns out that Eq. (A9) in the light-front coordinates 

results in a time-dependent equation only for two components among the four. For the other two components, 
it gives a time-independent constraint. The elimination of these two components from the Hamiltonian, by 
using this constraint, generates an extra term in the Hamiltonian, which is just the contact interaction. 

To demonstrate this in more detail, let us split the four-component spinor * into two two-component 
pieces, *(+) and defined as 

*(±) = A(±)*, (AlO) 
where A(='=) = ^7^7^ with 7^ =: 7° ± are projection operators with the properties 

A(±)2 = A(±), A(+)A(-) = A(-)A(+) = 0, A(±)^ = A(±), 



'*In principle, one can carry out the derivation without specifying the particular four-vector oj = (1,0,0,-1), but 
introducing four orthogonal four- vectors, one of which is ui. 
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The equation of motion (A9) also splits into two equations 

ja-^r(+) _ (ia-L-a-L + TO7° - ff$7°)*(-) = 0, (All) 

- {ia-^-d-^ + m7° - 5$7°)*(+) = 0. (A12) 

It is important that Eq. (A12) does not contain the "time" derivative d~ = "^-^^ and, hence, represents a 
constraint connecting the components and at any "time": 

= -^{ia^-d^ + m7° - 5$7°)*(+). (A13) 

The Hermitian integral operator acts on a function / according to the following rule: 

^/(^-) = -\j dy-e{x- - y-)f{y-) , (A14) 

where e{x) = 1 for x > and e{x) = — 1 for x < 0. In momentum space it corresponds to the division by 
(or by u)-p), where p is the four-momentum conjugated to the coordinate x. 

We represent the state vector in Fock space formed by the action of the usual (free) creation operators 
on the vacuum. One should also express the Hamiltonian in terms of these free operators, i.e., the operators 
in Schrodinger representation. The momentum operator P~ in Schrodinger representation has the same form 
as V~ , changing Heisenberg field operators to Schrodinger ones. The Heisenberg operators ^'("'"^ and $ are 
connected with the corresponding Schrodinger operators V^"*"^ and by the relations 

^(+) ^^gip-a.+ ^g-ip-x+^ (A15) 

where it is supposed that at a;+ = the Heisenberg and Schrodinger operators coincide. As it can be seen 
from Eq. (A15), the possibility to express the momentum operator through the free fields is essentially based 
on the existence of a "time" moment (chosen to be .x+ = 0), when all field components entering this operator 
are free simultaneously. However, in CLFD the constraint (A13) connecting the spinor field components \['(+) 
and ^f^") at any time moment involves interaction terms. By this reason in order to construct the momentum 
operator one must express the latter through independent field components only and after that substitute the 
Heisenberg field operators by the Schrodinger ones. In the case considered above the independent components 
are and 

So, changing in Eq. (A8) by ■^^"'"^ and $ by y, denoting 



id+ 



(+) 



where ^ satisfies the free Dirac equation, and returning to the covariant notations, we find for the interaction 
part of the four-momentum: 



j 6{uj-x)d^ 



P^J'^ = uj„ I 5(u;-x)(fx 



(A16) 



The integrand in Eq. (A16) is nothing else than the interaction CLFD Hamiltonian. 

To get the complete Hamiltonian incorporating the mass renormalization and the specific CLFD coun- 
terterm discussed in Sec. Ill [see Eq. (27)] we should add to Eq. (Al) the following additional terms: 

(5m** + Z^^ I ^ - 1 ) * EE + Z^^^^^, (A17) 

\ tu)-o J iLO-a 

where 5m' is defined by Eq. (33). It is easy to see that the counterterm with Z,^ in Eq. (A17) just generates 
Eq. (27) in momentum representation. 
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A question may arise whether it makes sense to add the counterterm with Z^^, expUcitly depending 
on the light-front position, to the Lagrangian. On the one hand, the Lagrangian as a fundamental quantity 
describing the physical system considered may "know" nothing about CLFD and truncated Fock space. On 
the other hand, we can not neglect the fact that the structure of counterterms substantially depends on the 
method of calculations. Since CLFD is based on Hamiltonian dynamics, and the Hamiltonian is derived from 
the Lagrangian, wc can incorporate necessary counterterms, formally including the lattcrs directly in the initial 
Lagrangian. At this level such a procedure looks like a formal trick allowing to simplify the construction of the 
Hamiltonian. 

After adding the counterterms (A17) we need not repeat the derivation of the Hamiltonian from the very 
beginning. In a first step, we put .^o; = on the r. h. s. of Eq. (A17) and take into account the term with 5m' 
only. The additional term 5m!^^ in the Lagrangian changes the equation of motion (A9) to 

(z ^-m)^' = -(c?$ + 5m')*. (A18) 

Note that Eq. (A18) can be obtained from Eq. (A9) by the simple substitution <i> ^ $ + 8m' j g. Hence, 
making an analogous substitution ^ + 5m' /g in Eq. (A16), we find the mass operator including mass 
renormalization effects. In a second step, we add the counterterm Zi^^[m t/j / i{uj-d)]^ to the Lagrangian. Since 

^ = 2\l/(+)^\l/*^+), this counterterm involves the \E''^+' -components only. By this reason, it does not affect 
the constraint (A13) and comes into the Hamiltonian as an addendum, without changing its initial form. 

Summarizing, we conclude that the final momentum operator including the whole set of counterterms can 
be obtained from Eq. (A16) by changing ^p ^ ip + 5m' / g and by adding to the resulting integrand the quantity 
—Z^-^lm i/)/i{ui-d)]tp. After simple algebraic transformations we get the regularized momentum operator 

p;"* =u;pj d^x5{u;x) W\x), (A19) 

with i?'"*(a;) given by Eq. (34). Except for the "normal" interaction —gijjip<p, all other terms in the Hamilto- 
nian (34) are contact ones. 

2. Vector boson 

The Lagrangian of a system of interacting spinor and massless vector fields is, also at the tree level, 
£ = 1 [^Yid.'^) - (5.*)7"*] - m** - ^Fi^^F^, _ ^(a . $)2 + (A20) 

where F'^'^ = 9^cf)'^ — 9''$^. Wc have added to the Lagrangian a term proportional to the square of the vector 
field divergence in order to have a possibility to choose different gauges. Due to the gauge invariance all physical 
results must be independent of the choice of A. 
The equations of motion read: 

(i ^ - m)* = -g (A21) 

d^^" - (1 - X)d''{d ■ $) = -5*7"^*, (A22) 
where d'^ = d ■ d = d^d~ — . The momentum operator acquires the form 
p- = ^ J (fx^dx- |-i*7-a+* + 2i*7^ • d^-^ + 2m** 
+1^ [(9-$+)2 - [d+^-f] + \{d^ ■ ^^f + [d^ X 

-{8^^-) ■ (a-^$+) + (1 - A)(9^$-) • (9+*-^) - 2g* . (A23) 
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Analogously to the scalar boson case, only two spinor components coming into are independent, \J/(~) 

being expressed through To demonstrate it we split, as before, the spinor ^ into and Instead 

of Eq. (A12) we now have 



{id+ + - (ia-L • + rn7° - 57°AW ^)^-(+) = o. 



(A24) 



Note that since A(+) ^^'(+) = -7-'- • the $~ component does not contribute to Eq. (A24) at all. The 

solution of this equation is 



(A25) 



(-) 



id+ 



(A26) 



10+ iV+ 



where the Hermitian operator is defined as 



iV+ 



fix-) = -'-j dy-e{x- - y-) exp 



f{y-)- 



(A27) 



(A28) 



The operator is constructed so that it turns to zero for any f{x~) if <? = or $+ = 0. 

As far as the vector field is concerned, the situation differs strongly from that for scalar bosons. Indeed, 
the vector field has not one, but four components satisfying Eq. (A22). In principle, some restrictions can be 
imposed on these components, reducing the number of independent ones (this just corresponds to a certain 
choice of gauge). Besides that, the momentum operator (A23) depends on the first order "time" derivative 
while Eq. (A8) does not depend on 9~<I>. When constructing the momentum operator in Schrodinger 
representation, we can not merely substitute d~^~^ by the corresponding free field derivative d~ip~^. The reason 
for this is that is not an independent quantity, because Eq. (A22) taken for i/ = + allows to express 

through the other field components and their "spatial" derivatives: 



1 + A 



2d^ ■ *-^) 



It is important that this connection includes interaction terms. 

Applying the operator di, to the both sides of Eq. (A22) gives 

A52(S- $) = 0. 



(A29) 



(A30) 



In deriving Eq. (A30) we took into account the current conservation, 9,^(^7'^'!') = 0, which, in its turn, follows 
from Eq. (A21). In spite of that, Eq. (A30) contains vector field components only and cannot be regarded as a 
constraint because it involves "time" derivatives of the second order. 



a. Feynman gauge 



In the Feynman gauge all the four vector field components are treated as being independent. The light- 
front free vector boson propagator is — g ^^5{q^)6(u!q) . To meet these requirements one should put A = 1. Indeed, 
consider the solution of Eq. (A22) corresponding to a free boson (i. e. at g = 0). This solution is a superposition 
of plane waves (p''(g)e~''^ with = 0. Substituting it into Eq. (A22), one gets (1 - X)q''[q ■ ip{q)] = 0. If 
Q ■ fil) = 0; only three field components are independent, but not four. So, we are forced to put A = 1. Under 
this condition, the momentum operator (A23) does not depend on d~^~^. 
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Since in the Feynman gauge the independent fields are and the subsequent analysis is quite 
similar to that made for the scalar boson case. It leads to the following expression for the interaction part of 
the momentum operator in Schrodinger representation: 



jint 



with the interaction Hamiltonian 



2i{oj ■ d) 



■9 ^ 



2i{u ■ D) 



(A31) 



(A32) 



The operator ■ D) is given by the r. h. s. of Eq. (A28) after the replacement 
explicitly covariant form it can be rewritten as 



1 



i{uj-D) 



f 



R* 



1 



-R 



i{uj-d) i{uj-d) 



f, 



(f^ = ui ■ (f. In an 



(A33) 



where R = exp • f)] ■ Note that the operator jfjjrg^ inside the exponent acts on u!-(p only, while the other 

such operators in Eq. (A33) act on all functions standing to the right of them. 

The second term in Eq. (A32) is a fermion contact term. In the Feynman gauge, to order g'^, contact 
terms related to the massless vector field are absent. The decomposition of the last term on the r. h. s. of 
Eq. (A32) starts with and contains an infinite number of terms. The matrix element of a term of arbitrary 
order contains the contractions of the vector field operators entering Eq. (A32) between themselves and with the 
external fields. As is shown in Sec. VA, in the Feynman gauge the internal contractions disappear since they 
are always proportional to = 0. The contractions with external fields only survive. Therefore, in truncated 
Fock space, the number of contact terms is finite and the Hamiltonian (A32) has polynomial dependence on the 
coupling constant and field operators. We emphasize that this property is a merit of the Feynman gauge, while 
another choice of gauge may lead to appearance of an infinite number of contact terms even for truncated Fock 
space. 

Incorporating the self-interaction Sm'^"^ by means of the substitution ^ ^^ + 6m' /g [with 5m' given by 
Eq. (33)] and adding the counterterm —Zij'tplm f/[)/i{u}-d)]ip, we obtain the Hamiltonian (71) with H'{x) given 
by Eqs. (72)-(73c). 



b. Light-cone gauge 

Maximal simplifications arc obtained in the LC gauge where — (equivalent to ti; • $ = in CLFD). 
In this gauge the bosonic propagator —[g^v — {qui^v + QL/'^/i)/(w • q)]5{q'^)9{u! ■ q) is transversal to to in both 
indices. Note that expressions for vector meson propagators in the LC gauge can be found in Refs. [2,3]. 

For $+ = 0, we get from Eq. (A29): 



4,9 



1 - A 



J_*(+)^^(+). 



(A34) 



This is a constraint indicating that only two transversal components among four initial ones, , are inde- 
pendent. Under these requirements the vector field, generally speaking, does not satisfy the Lorentz condition: 



23 



L^(+)^'$(+) ^ 0. 



The relation (A30) now reads: 



^2 



1- A 



d+ 



= 0. 



(A35) 



(A36) 



To meet the last equation one has to put A = 0. 
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Now we can take Eq. (A23) and set A = 0, = 0. Then, using Eq. (A34) we get the momentum 
operator expressed through the independent fields In Schrodinger representation we obtain 



5mt 



(A37) 







2i{u-d) 



1 



(A38) 



The first term in Eq. (A38) is the standard spinor-vector interaction Hamiltonian, while the second and the 
third ones are the contact terms. Adding to this Hamiltonian all necessary counterterms, one gets Eq. (71) with 
H'{x) given by Eq. (93). 



APPENDIX B: CALCULATION OF COUNTERTERMS 



1. Scalar boson 



Explicit calculation of the fermion self-energy (15) is similar to the one given in Ref. [8]. We will also 
calculate the contribution (16) of the fermion contact term. From Eq. (15) we can express the functions A{p^), 
B{p^), and Cip^) through S(p): 



g^Aip") = iTr[E(p)], g^Bip-") = -^Tr[S(p) 0], g^Cip") = -^-Tr 
4 4a; -p 4m 



E(p) 



uj-p 



(Bl) 



For our purposes it is enough to know these function in the region > only. 

We introduce the variables fc+ = fco + ^z, = {kx^ ky) and x = k'^ jp'^ . For convenience, we will carry 
out calculations in the frame where p = and uj = (1, 0, 0,-1). Hence, 



k+ 



,2' "'-^ ~ 2k+ ' °~ ''2" 



Substituting the integral (15) for E(p) in Eq. (Bl) [and similarly for Eq. (16)], we get 



A{p') 
Bip') 

Cfc 



m 
m 



dkl 



dx- 



1 



k'j_ + m^x — p'^x{l — x) + /u2(l — x) ' 

/"^ l-x 
Jo k]_ + m'^x - p'^x{l - x) + ^'^{1 - x)' 

rl-S 



k]_ +vnf -p^{l - xf 



■327r2myo {\- x)\k\+m^x-p■ 

A_ f dk' \ dx ^ r ^ 

2'K'^m Jo Js x{l - x) Ji+s x{l - 



■x{l — x) + /x^(l — x) 



(B2a) 
(B2b) 
(B2c) 
(B2d) 



We introduced in Eqs. (B2a)-(B2d) two cutoffs L and S in the variables kj_ and x, respectively, which restrict 
the integration region as 

S <x <l-6, l + ^<x<+oo, fci < i^. 

In the following it is implied that L 00 and i5 ^ 0. The two cutoffs do not constrain the integration in a 
spherically symmcitric domain and, in principle, violate rotational invariance. We will see that 6m depends on 
the cutoff L logarithmically and does not depend on 5. 

In order to calculate the constants A, B, and C we use Eqs. (B2a)-(B2c) at p"^ = rv?. We thus find 
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m f^^ 1 
A = --—r I dk\ / dx—, — -, (B3a) 



771 


r 


167r2 


I 


m 




167r2 


I 


1 





B=-T^,i dkii (B3b) 

- ^, 2 ^ + m''x{2 - x) .^^ , 

C = / / dx— ttA ;r-^ ijz TT. (B3c) 

Calculating the integrals (B3a)-(B3c) for A, B, C, and those for C/c in Eq. (B2d), and retaining only the terms 
divergent at i — > oo and ^ ^ 0, we get 

^ = -T^log^, (B4a) 

B = -^,^og-„ (B4b) 



1 R 2 


^2 


m , 
32.^ 




X2 log 5 


m 


327r2m 


327r' 


L^log^ 




327r2m ' 





C= ^:r^-:j^log^, (B4c) 



Prom here we obtain 



and 



2 r 2 

5m = g''{A + B + C + Cf,) = -^-^\og—. (B6) 

Note that the strongest singularities [oc a^fj L^log^] cancel in the sum C + C/c- 
The normalization integral (65) reads 

2 2 f 9^ Z"^' „2 /"^ f'fci +m2(,x-2)2la; 1 



The leading term of A'' at L ^ oo is 



^=^log-^. (B8) 



Now let us prove that C{p^) given by Eq. (B2c) is a constant. For this purpose we calculate the difference 

AC(p2) = C(p2) - C(m2). (B9) 

Substituting C(p^) from Eq. (B2c) into Eq. (B9), we obtain 

2^ _ffV^m^ /• .2 (' [kl{2x-l) + m'x^-y.\l-xf]dx 

' 327r2m 7 ^ [kl + - (1 - x){p'x - iJ?)][kl + m2x2 x)] " ^ ^ 

The integrand has no singularities at a; = or a: = 1. Integrating over x results in zero for arbitrary p^: 

ACip'^) = 0. (Bll) 

This means that the w-dependent part, Eq. (32), of the fermion-boson scattering amplitude disappears, and the 
whole amplitude, Eq. (20), becomes w-independent. 
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2. Gauge boson in the Feynman gauge 



The fermion self-energy with the gauge boson propagator in the Feynman gauge is given by Eq. (78). 

From this equation follows that the quantity S^j, differs from the on-mass-shell self-energy for the scalar boson 
case, S, defined by Eqs. (52) and (15), by the factor —g^i, and by that the gauge boson mass is zero. Hence, 
the coefficients Ai, Bi, and Ci in Eq. (79) can be obtained directly from Eqs. (B4a) (B4c) by setting /i = 0, 
that does not make any influence on their asymptotics' at L ^ oo. Then, using the relations (80), we can easily 
find the coefficients A, B, and C determining the self-energy (78) in the Feynman gauge. Therefore without 
repeating our previous calculations we obtain 

C= ^-T^^log^, (B12C) 



This gives 



^.=fl^(C + C,J = -01og^, (B13) 



Sm = g\A + B + C + Cf,) = ^log^. 



(B14) 



Like in the scalar boson case, the terms quadratically divergent in and depending on the cutoff 6 cancel in 
the sum C + Cfc- 



3. Gauge boson in the Light-Cone gauge 

The self-energy in the LC gauge is given by Eq. (98). Doing the same as in deriving Eqs. (B3a)-(B3c), 
we obtain for the coefficients A, B, C [and similarly for the contributions from the fermion and boson contact 
terms Cfc and C^c defined by Eqs. (94) and (96), respectively]: 



8^/ *i 







dx 



1 



B = - 



m 
8^ 



dki 



kj_+ rrfix'^ 



dx 



kj + mP'x'^ 



C 



Cfc = 



Cbc — 



dx 



k^ipc^ -2x + 2)+ m'^x^{2 - x) 



x^{l — x){k\ + wp-x^) 
IQ-K^m Jo Js x{l - x) Ji+s xil - x) 



rl-5 



r pi-s 

dkl 

Jo Jo 



dx 



{1-xy 



dx 



{1-xy 



dx 



{1 + xy 



(B15a) 
(B15b) 
(B15c) 
(B15d) 
(B15e) 



Calculating these integrals, taking the limits L — > oo and 5 — > and keeping divergent terms only, we get 
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Prom here we find 



and 



Sw-'m loTT^m STr^m Idtt^ 
Qc = -^, (B16d) 

SiT-'m STT^m 



= + + Cfce) = log ^ (B17) 



2 r 2 

= ^^(^ + B + C + Cfc + Cbc) = ^ log ^. (B18) 

The terms C, C/c, and Cbc are quadratically and logarithmicaUy divergent in the cutoff L and divergent in 
the cutoff S like 1/6 and log^. However, aU the divergences excepting log(L^/m^) ones cancel in the sum 
C + Cfc + Cbc- 

APPENDIX C: EQUATIONS FOR THE SPIN COMPONENTS 

1. Scalar boson 

The spin components ai, bi, and 62 in the scalar boson case are defined by Eqs. (53) and (56). As 
explained in Sec. IV C, the equations for the spin components are obtained by substituting these equations into 
Eqs. (54) and (55) and taking traces. As a result, we get the following system of three homogeneous equations 
for the three unknown constants ai, 61, and b2'- 

2aim{5m'^ - 2m6m' - 2mZ) - big[{A + B){Sm' - 2m) - 2mC] 

+ 2b2g[{A + B)m - 6m' B] = 0, (CI) 

2ai5m[m^(3 + x)- m5m'{l + x) + s(l - x)] + bi[g'^m{A + B){1 + x) - m^(3 + x) - s{l - x)] 

+ 2b2m{g^B -m){l+x) = 0, (C2) 

2aigm[m{l + x) - dm'x] + bi[g^{A + B)x - m(l + x)] + 262(3^^ - 'm)x = 0. (C3) 
The condition that the determinant (59) of this system is zero results in the solution (60), (61). 



2. Feynman gauge 



In the case of gauge boson in the Feynman gauge, the one-body fermion vertex function Fi is determined 
by one spin component only, as given by Eq. (84), while the two-body vertex function F2 by eight spin compo- 
nents, four of them being zeroes when sectors higher than \Fb) are omitted in calculations. The derivation of 
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the system of equations for the spin components is described in Sec. VC. The system of five equations for m, 
bi-4 is given by 



aim{2mdm' - 6m' + 2mZ) - big[{2Ai - Bi){Sm' - 2m) + 2Cim] 

+ b2mg{Ai + Bi) - 2hzg[Aim + Bi{5m' - 2m)] + 2bigBim = 0, (C4) 

2a\gm [m^(3 — x) — s(l — x) — mdm.'{2 — x)] 

- bi [2mg^{2Ai - Bi){2 - x) + {?, - x) - s{l - x)] 

- b2m?{l + x)+ 2bzm[2g^Bi{x - 2) + m(l - 2x)] - 264^^^ = 0, (C5) 



{2aigm - bi){l + x) = 0, 



(C6) 



2aigm[x5m' + m(l - 2a;)] + bi[2g'^{2Ai - Bi)x - m(l - 2x)] + 2b3{2g'^Bi - m)x = 0, (C7) 



2aig {m^(l +x) + \2m^ + s(l - x)\(m - 5m')} 



6Jm^(l + a;) + 



[2to2 + s(l - a;)][2£f2(2Ai - Bi) +m] 



m 



bi\m^(2, + x)^ s{\- x)\-2bi 



n?{\^x)+2g^Bx 



2rv? + s(l - x) 



m 



2b4m^{l + x) = 0. (C8) 



The determinant of this system of equations is given by Eq. (87). The condition Det = results in the 
solution (88)-(90). 

3. Light-Cone gauge 

The equations for the spin components, derived in Sec. VIC, have the form 

aim[2mSm' - 6m^ + 2mZ] + big[5m'{Bi - Ai) + 2m{Ai - 2Bi -B2 + 5Ci + 'iC[)] 

+ 2b2g[m{Bi - Ai) - BiSm'] = 0, (C9) 



2aigm 



1 + x' 



1-x 



(s — m^) — 2m^ + m6m'{l — x) 



-61 



1 + x' 
1-x 



(s - m^) - 2m^ - 2g^m{Ai - Bi)(l - x) 



2b2m{2g^Bj_-m){l-x) = 0, (CIO) 



2aigm[m{l - x) + x6m'] - 6i[m(l - x) - 2g'^{Ai - Bi)x] + 2b2{2g'^Bi - m)x = 0. (Cll) 
The condition that the determinant of this system is zero [see Eq. (106)] results in the solution (107)-(109). 



35 



APPENDIX D: DEPENDENCE OF 5M ON THE TYPE OF CUTOFF 



Here we will show how to adjust the cutoffs in Eqs. (B15a)-(B15e) in order to reproduce the perturbative 
value of 5m obtained in Ref. [6] for the LC gauge. 

The coefficients B, and C proceed from the light-front diagram shown in Fig. 2(a), whereas Cfc and 
Cbc originate from the loop diagrams shown in Figs. 2(b) and 12, respectively. Following Ref. [6], we regularize 
the amplitude of the diagram of Fig. 2(a) by the condition 



s = — + 

X 1 — X 



with the additional requirement 

fc_L < i < A. 

Eq. (Dl) leads to the restrictions a < x < 1 — l3, where 



(Dl) 



(D2) 



(D3) 



The integrands in the expressions for A and B [see Eqs. (B15a), (B15b)] are not singular in x, hence, A and B 
are given by the same formulas (B16a), (B16b). The expression (B15c) transforms to 



C 



1 



IGTT^m Jq 



^ „2 f^~\ kW -2x + 2) + m^x'^{2^x) 
dkj_ I dx - ' 



1 



dkl 



a;2(l — x){k']_ + mPx^) 

2m?{l - x)dx 2dx f'^'^ dx 

Jq m?x'^ + fc^ Ja Jo 1 - ^ 



m ^ L _^ 1 



dkl 



log/3-- + 2 
a 



(D4) 



As far as the divergent integrals over x for the coefficients Cfc and Cbc in Eqs. (B15d), (B15e) arc concerned, 
the prescription of Ref. [6] demands to use for their regularization the only cutoff a. Therefore, we have 



Cfc = 



IGtt^to Jo 



dkl 



dx 



x{l — x) 



dx 



1+ 



a X{l-X)^ 



1 



dkj_ log a, 



(D5) 



- 1 ^1.2 \ f^~" dx f°° dx _ f 

"-'"-le^^mjQ *^Uo {i-xr+J,^c.a-xr Jo 



dx 



[i + xy 



1 



IQiT^m Jq 



[ dkl 




Jo 


a 



Now, using the first of Eqs. (Ill) we get 



Sm ■ 



■ log + 



[ dkllog^. 

Jo OL 



Taking into account Eqs. (D3) and performing the integration we finally obtain 



(D6) 



(D7) 



(D8) 
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The same coefficient 3/16 as in Eq. (D8) appears also for the regularized perturbative SrriiDF in standard 4DF 
approach: 



5m4z,^ = ^log%, (D9) 



where L^df is a relativistically invariant cutoff restricting the region of integration in four- dimensional space. 
Note however that in spite of the similarity between Eqs. (D8) and (D9) the cutoffs L and Lajp restrict different 
integration regions and therefore have different meaning. 

So, the only difference between our method of calculation and the method from Ref. [6] consists in that 
integrating over x we apply the same cutoff 5 for all the divergent integrals in Eqs. (B15c)-(B15c), while the 
authors of Ref. [6] used two different cutoffs, a and /3 (dependent on k\), as exposed above. Putting a = /3 in 
Eq. (D7) immediately reproduces our result (B18) for 5m. 
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FIG. 1. Fermion-boson CLFD scattering amplitude in p^-order of perturbation theory 
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FIG. 2. Mass operator in CLFD. The lowest order contributions are shown: the light-front self-energy, — (a), and 

the contact term, — E/c(p) (b). Here p = pi — wn. 
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FIG. 3. Interaction vertices generated by the Hamiltonian Hi{x), Eq. (35a). 
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FIG. 4. Interaction vertices generated by the Hamiltonian Hic{x), Eq. (35b). 




FIG. 5. Interaction vertex generated by the Hamiltonian H2{x), Eq. (35c). 
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FIG. 6. Interaction vertices generated by the Hamiltonian H2c{x), Eq. (35d). 
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FIG. 7. Interaction vertex generated by the Hamiltonian H2c{x), Eq. (35e) 



FIG. 8. Interaction vertex generated by the Hamiltonian H3{x), Eq. (35f). 
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FIG. 9. System of equations for the one- and two-body vertex functions for the scalar boson case. 
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FIG. 11. System of equations for the one- and two-body vertex functions for the case of vector boson in the Feynman 
gauge. 
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FIG. 12. Interaction vertices generated in the LC gauge by the piece H'j^c{x), Eq. (93), of the light-front Hamiltonian. 
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